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ABSTRACT 


This research 185 a study in seepage analysis for engineered 
structures built below the groundwater tables in particular an 


analysis of Underdrain systems for such structures. Ganeral 


seepage theories and design cansiderations are discussed as 
part of the literature review. Alsa design concepts and 
methods of analysis are covered which gives a brief overview of 
design considerations and techniques which are currently used 


for this type of seepage problem. 


Rey variables to the design of an underdrain system are stuclied 
im detail by using a two dimensional, steady state finite 
element model. Such variables studied are drain spacing, drain 
width, anisotropic hydraulic conductivity, thickness ot the 


aquifer, head differential and sloping aquifers. 


Results are presented in the farm of dimesionless plots and 
noamographs for drain discharge and free surface plots for the 


different variables studied. 
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CHAFTER ONE 


INTRODUCTION 


For civil engineered structures which interact with soil, a 
major consideration for design is the presence of oroundwater. 
Dams, buildings, dry dock facilities, underground structures, 
pavements and earth retaining walls must be designed 
considering the quantity of seepage, the position of the 
free surface, seepage forces and seepage control. Good 
drainage incorporated in the design of these and other 
structures can often benefit through substantial cost 
reductions as well as increasing the safety and usefulness of 


the structure. 


Many embarassing and sometimes disastrous problems oaften result 

from man’s endeavors to control seepade. Saturated soil 

backfills have caused many retaining wall failures aque to 
increased hydrostatic pressures following extensive rainfall. 
Highways and roadways prematurely deteriorate and fail due to 
rapped water in the pavement section. Basement slabs, pond 
liners and drydock facilities Nave been severely damaged by 
Hydrostatic uplift pressures and seepage. Most dam and 

Slope failures can be attributed to excessive seepage or 

uplift pressures. These examples and untold others indicate 


that good seepage and drainage control together with an adequate 





understanding of groundwater 1s a key element in the design of 


most civil engineering structures. 


Many Civil engineering structures are located in areas where 
groundwater causes problems of uplift and seepage. The design 
therfore must consider all problems associated with high 


~ 


Br aupawat er and what measures must be-taken to assure the 
integrity of the overall design. Casagrande(1965), in talking 
about engineering risks, concludes that, "The margin of safety 
that we incorporate into our structures should bear a direct 
relationship to the magnitude of potential losses, and it must 
also take into account the range of uncertainty involved." 


This important concept hopefully will be considered throughout 


the presentation of facts and data obtained in this paper. 


Frior to the 20th century, analysis of seepage and the control 
Of drainage for enginerring structures was done by experienced 
engineers who based their designs on trial and error methods 
which sometimes were successful and sometimes not. Though 
experience has proven to be very neccesary and helpful in the 
design of drainage systems, very useful tools employed today 
are numerical methods such as finite element and finite 
difference techniques which require the use of the digital 
computer. Their problem solving capabilities have made 
possible solving very complicated and complex groundwater 


problems. 





The purpose of this research is to study the design of 
underdrain systems for large structures placed below the 
groundwater table. In chapter two, basic seepage theories and 
design requirements, designs concepts and useful methods of 
analysis are presented as a Brief literature review. This 
chapter Pavers critical aspects of designs that should be 


considered and contains a brief review of methods used by 


geotechnical engineers. 


In chapter three, the finite element program developed by Dr. 
M. Aral, of the School of Civil Engineering at Georgia 
Institute of Technology, 1s discussed describing how the 
design of Building underdrain systems can be modeled using 
this technique. This finite element program made possible 
studying effects of key variables such as permeability, 
hydraulic head, drain geometry, depth of groundwater, and 
anisotropic soils. Frogram verification is performed using 
another program, SEEFS4, modified by the Civil Engineering 
Department at Virginia Folytechnic Institute and State 
University, for the microcomputer. Dupuit*’s solution is also 
used to compare selected cases with the finite element 


programs. 


fy 





Chapter four presents data generated from the finite elament 
study such as position of the free surface and flow quantity 
for selected systems of underdrains suitable for use beneath a 
foundation. The data is presented in the form of charts ard 


nomographs to be used as a guide for underdrain system desiqn. 


Chapter five summarizes the results of the study and presants 


appropriate conclussions. 


Appendix A contains program documentation for the finite 
element program used in this study which details to the user, 
documentation to the use of the program. Appendix B presents 


example input and output data generated by this study. 





CHAFTER 2 


LITERATURE REVIEW 


eel Seepage Theories and Design Requirements 
\ 


Sroundwater 


SOL1S, sametimes described as a porous media, 15 a mass made up 
of solid matter having interconnected voids, passages, fissures 
or pore spaces. These pore spaces can be filled by a liquid 
and or a gaseous phase. Figure 2.1 shows several typical soil 
and rack indices which depicts how pore spaces can act as 
conduits for water or subcapillary openings in which water may 


be held by adhesive capillary forces (Me@inzer,1942). 


The soail mass forms a complex network of irregular passages 
through which the fluid may flow. A key parameter of 
groundwater movement is the porosity of the soil. AS gravity 
acts on the fluid, it moves through the pore spaces at varying 
rates through each apen pore space. From a macro scale, this 
flow has been mathematically described by several authors 
(Bear, 19681; Terzaghi, 19423; Cedergren, 1977s Harr, 1962; 


Bouwer, 1973;Freeze, 1971). 
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Figure 2.1 Diagram showing soil and rock interstices. A. Well 
sorted sand deposit having high porosity: &. roorly sorted 
deposit having low porosity: C. Farticles that may be porous 
themselves; D. Deposit whose porosity may be diminished by 
deposition of mineral matter in the interstices; E. Rock that 
Remenols by Solutions F. Rock that is porous by fracturing 
(after Meinzer, 1942) 
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Moisture content in soils can vary with depth from the surtace, 


as shown by Figure 2.2. Most seepage problems considers the 


= 


iL 


of saturation which is the zone where the soil mass 15 
completely saturated. The terms groundwater level, phreatic 
water Level, water table, free surface or groundwater surface 
are all synonymous to denote the level of water where 


atmospheric pressure exists. 
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Figure 2.2 Subsurface groundwater distribution. 
(after Bear, 1979) 


In dealing with capillary fringe in fine grained soils the 
height of water can very significant. As the pore diameter 
decreases, the capillary rise increases. Although the soil is 
considered unsaturated, the soils moisture comtent will be 


higher with depth. 


} 


The zone of capillary rise and fringe zone may contain large 


quantities of water for fine grained soils. Water removal by 
drainage, due to capillary forces, may be difficult for such 
fine grained soils. Figure 2.232 shows how the actual moisture 
distribution is approximated by a step distribution. This 
approximation of the capillary fringe level above the free 
surface, Ne, 1s Often neglected in some seepage problems. 
Several emperical methods, based an orain size diameter and 
porasity of soll, are available to estimate ae tad eet) 


foeiemet7o7. Folubarinova, 1952, 1962s Silin Bekchurin, 1958). 
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fringe zone. (after Bear, 1979) 


Figure 2.3 Approximation of phreatic surtace and capillary 


Many textbooks on seepage and groundwater discuss water bearing 
and retention capibilities of the capillary fringe and 

movement of moisture through this zone (Terzaghi, i942; EFovacs, 
feos ouWwers, 1979). However, when this height of water 185 very 
small relative to the total thickness af the aquifer, the 
capillary fringe 15 sometimes neglected. The geotechnical 
engineer should be concerned with the unsaturated-saturated 

flow conditian when the stability of the engineered structure 


is to be analyzed, such as a dam. 


Groundwater analysis for seepage problems must consider the 
quantity of water to be drained, the soils resistance to 
drainage, and what effects the drainage may have on the soil 


and the structure. Movement of groundwater can be categorized 





either as seepage; capillary rise or tlow; percolation; or 
turbulent flow. For this study, only movement of groundwater 


by seepage is considered. 


Design Requirements 


Structures built below the watertable must be designed 
considering the following: 

(1) structural ability to withstand the applied hydrostatic 
water pressures, 

(2) drainage control which short-cuts the flow or movement 
of groundwater from getting to the structure, 

(2) drainage control which is incorporated with the overall 
design to remove the groundwater, 

(4) assure water velocities be small enough to prevent 


excess particle migration and erosion. 


Savings of dollars, time and natural resources can be 
recognized when structures can be considered on dry, stable 
materialsin addition to thinner basement slab thicknesses may 


be used. 


10 





Darcy’ s Law 


In 1856, Henry Darcy investigated the flow of water in vertical 
homogeneous sand. This work has developed the most fundamental 
relationship which is still used in seepage analysis today. 


Darcy*s Law states: 


where @ (L3/T) is a volume discharge or flow of a given 
cross-sectional area, A en in a given time. The discharge is 
also a function of the hydraulic gradient, i (L/L) and a 
resistance coefficient, k (L/T), the permeability of the soil. 
Investigators, Muskat (1927), Taylor (1949) and Leonards 
(1962), have concluded that Darcy’s Law is valid where 


velocities are low or when flow is laminar, i.e. Reynolds 


number less than 190. 


Hydraulic Conductivity 


The soils resistance to flow or the coefficient of 
porportionality, k, for Darcy’s Law, 15 also called the 
hydraulic conductivity. Hydraulic conductivity is a 
coefficient which expresses the ease Or resistance the water 
has flowing through the soil. It depends on the soils matrix 


and fluid properties of the permeant. 


IME 





Many researchers have performed tests to attempt to develope 
empirical relationships of the soils properties to hydraulic 
canductivity. Some of these properties considered are grain or 
pore size distribution, shape of the grains or pores, 

tortuosity, specific surface and porosity (Mitchell, 197463; Eozeny, 
1927: Carman, 1956; Cedergren, 19773; EFovacs, 1901). All report 
that careful field investiagation and careful conductivity 
etudies are necessary for each design. Figure 2.4 is presented 
to show some characteristics of variaqus soil types and most 
important, methods for determination of selected properties 


(Terzaghi and Feck, 1°48). 


In the unmsaturated zone, the hydraulic conductivity, moisture 
content and pressure head is extremely complex and not well 
understood. Again, field measurement is the best method for 
abtaining the best parameters, such as the use of theistor 
psychrometer (Fay and Low, 1979) or gamma ray absorption 
techniques (Ferguson, 19703; Aral, 1963). However, many 
empirical formulas have been developed to relate the soil 
parameters to a relative hydraulic conductivity of the 
un-saturated zone (Kovacs, 1981; Irmay, 19534; Childs and 
Collis-George, 1950; Gardner, 1958; Brooks and Corey, 1966; 


Maslia, 1980; Reeves and Duriguid, 1975). 
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2. Design Concepts 
Seepage Forces 


The effects of water on the soil and the structure is 
controlled by the pressures in the water. In soil mechanics, 
the water pressure is defined as neutral stress. For water not 
mMOving, or steady state, the pressure can be computed from the 


basic equation of hydrostatics 


Foo+ g ne) 
Sw 


where h, energy head, is equal to the position, =, plus the 
pressure head or piezometric head, ie eee Typically in 
seepage designs, the velocity head is very small and may be 


neglected. Then the pressure head plus elevation is equal to 


head, h (Cedergren, 1977). 


For the condition that water pressure is greater than the total 
load, say of the soil and structure, the pressures are termed 
uplift pressures and liquifaction can occur. Therfore the 
control of uplift pressures must be a primary design 


consideration. 
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Men tical Velocity 


According to Darcy's Law, the discharge velocity of seepage 


through soil is 


seepage is only occuring through the pores and voids. 


Therefore the true seepage velocity should be a function of the 





poroesity of the soil mass. Re-writing equation (2.3) the 
seepage velocity can be expressed as 
k x i 
van = (2, 4) 
A 
e 


eee he ts the ertective porosity. 


In flow conditions, the saturation of soil 18 moving into or 
oot OM ehe Uneaturated soll. thejcapillary head, NM, ; 18 
working with gravity to increase the hydraulic gradient. For 


this condition the hydraulic gradient can be expressed 


— 
tJ 
Ch 


Figure 2.5 depicts this condition for flow through a sand. 


i 


rvolr 
Reservo rs 


iL 
Be oe a eo Us] | Reservoir 
Screen 
=(h+ h. iL 


Reference line 





meee «.2 Force of saturation. 


Roiling occurs when seepage pressures in an upward direction 
exceeds the downward force of the soil. This must be analysed 


ctl 


ri 


for all seepage designs. Critical hydraulic gradient, Le ; 


be expressed 


os 
- 
Hi 


where total unit weight of soil 


os 
tf 


UuNLt weight of water 


buaeyant unit weight of soil (NAVFAC,DM-7.1). 


os 
oO 
if 
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Filters 


Avery important design consideration is the possibility of 
seepage forces moving erodible sails and rock, termed piping. 
Betram (1940), under Terzaghi and Casagrande, established the 


following criteria for filter design, based on grain size. 


D (or filter) D (of filter) 
= C4tos <¢ 1s 


(of soil) o er 
D 95 of soil) Dis (of soil) ( ) 


This relationship is represented by Figure 2.6 to illustrate 


the prevention of piping by filters (Cedergren, 1977). 












LEGEND 
eer an. Dg; (soil) 
in-place soll a a Wi 
aor 
® = Dx; soil particle, cee 
entrapped tn filter NI 


[| = soil which has 
migrated into filter 
and is held by Ds; 
size soil particles 





Filter 
Nominal boundary 


before stabilization : 
under seepage eS 


eee. Lhhstrpation OF prevention of piping by filters. 
(after Rankilor, 1981) 
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reer aero Gf 045 Gf a filter to the Dge of a soll is 
considered the piping ratio. Many authors and organizations 
Pewee estanplished other criteria for filters, similar to that of 
Seuatvon (2.7) (355. Amemy Woerps of Engineers, 1941, 19535: U.S 


Bureau of Reclamation, 1°53; Sherard,etal, 1963; Cedergren, 


1977s Thanikachalam and Sakthivadivel, 19774). 


Filters must be designed to meet two major requirements; (1) 
they must be safe with respect to erosion and piping, and (2) 
must have sufficient discharge capibillities to remove the 
seepage quickly and efficiently. This must be done without 
inducing high seepage forces or Hydrostatic pressures. The 
right side of equation (2.7) states that the Pee Chas 
oaecieeosnculed fe At least + €O 4S trmes the Dj. of the soil. 
This insures that the hydraulic conductivity of the filter and 
drains 185 sufficient to prevent the Build-up of large seepage 


forces and hydrostatic pressures. 


Tnere is a wide range of engineering properties of materials to 
be considered for filter drainage systems desigqns. Froperties 
which may influence hydraulic conductivity are grain slze, 
density, mineralogical composition, the nature of the permeant oar 
the degree of saturation. Much work has gone toward methods 

for determination and factors which influence hydraulic 
conductivity (Milligan, 19753; Lambeand Whitman, 1969; 

Yemington, 1963; Houwer, 1978). There is still much research 


necessary to understand better this property of soil. 
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Trench Drains 


For the application of drainage for structures built below the 
Water table, the use of trench drains, which acts as collectors 
of the seepage or may contain collector pipes, 15 A very 
efficient method. Figure 2.7 shows the typical cross-section 
which may be considered for the application of this type of 
drainage system. Trench drains will be discussed further in 
the report and their utilization for this application will be 


studied in more detail. 


a BUILDING Ori 


FREE SURFACE —. + She 
TRENCH DRAINS 


Figure 2.7 Typical application of a trench drain system. 
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Geosynthetics 


Today, new synthetic membranes are used to serve the purpose of 


Gither filtration or seperation. 


important parameters for most seepage problems. 


This two factors are very 


Pigivie: 2.6 


shows comparisons of some of the features of membrane-wrapped 


drains to conventional granular pipe drains. 


ao 


Also Figure &. 


shows the bridging effect caused by good seperation of the 


larger materials and smaller particles. 


Granular drain 







Pipe drains 

need frequent 
rodding points. 
Expensive. 


Sand and stone must be 
supplied and laid to a strict 
specification. Very expensive. 


Pipe is expensive to buy, and 


Pipe must be must be laid by hand. Trench 


levelled by must therefore be wide, and 
engineer, to timbered if deep. 
be an exact 


gradient. This is expen- 
sive. Trench must be left 
open for this, and can get 


silted up during storms. unlike a pipe 


Well-point Wide trench 
pumping means more 

may be excavation and 
necessary to more backfill than is 
lay the pipe. technically necessary 


Figure 2.8 
membrane-wrapped trench drain. 
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The deep transporter means 
that settlement can take place 
without the drain blocking — 


Membrane-wrapped drain 






Cost of shoring 
up open trench 
is eliminated 


® 

= Granular fill only 
a Supports trench, 
3S so can be cheap 
ss specification 


This drain is very 
much quicker to 
construct & can 
be laid in flooded 
ground. 


a 


Can be 

taken deeper 
because sides 
do not need 
support 


Narrower 

width because no pipe to be 
laid, So men do not have to 
go into the trench. 


Conventional granular pipe trench drain and 
(after Rankilor, 1981) 





5 


Figure 2. 
material. 


Magnified cross section of permeable membrane 
filter drain 


drain 


aS 8; AD 
cbse. < 
ey 


bd 


ww? 
Py fod = nS 





Original Filter Bridging Drain 
soll zone network 
Structure in soil of larger 

(sometimes particles Membrane 

called 


Filter Cake) 





9 Bridging effects caused by use of 
(after Rankilor, 1981) 
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Stone filled 








membrane filter 





& good representation of how geosynthetics combine action of 


seperation and filtration, is shown by Figure 2.19. As water 


has more area of flow with the membrane, the seepage can be 


considered to be at a faster rate or drain quicker (Rankilor, 


heel ) 


Without membrane With membrane 


/ = i 
Water tries to 
escape through SECTION escapes 


small area rapidly 


Without membrane With membrane 


Small surface area Large stone/clay Large surface area 
for moisture loss contact area for moisture loss Stone/clay 
contact area 





Figure 2.10 Seperation effects By using a membrane. 
(after Rankilor, 1981) 
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a.o Mathods of Analysis 


In analyzing seepage problems; geometry, boundary conditions, 
Soil parameters, flow quantities and location of the free 
surface must be quantified By some means. Though this involves 
Many considerations, the engineer must gain reasonable safe 


results for the design. 


Through the history of seepage analysis many techniques and 

methods have been developed, such as graphical, analytical and 
numerical solutions. Still used today is the graphical method 
Of flow net construction. However with the complex geometries 
envolved and via the speed of computers, numerical models have 


been developed as very useful tools for the engineer. 


Flow Nets 


Many researchers have developed methods for analyzing flow 


through porous media. The oldest, developed by Casagrande and 


Terzaghi (1927), 18s the use of flow nets, a graphical approach 


to seepage analysis. 


Ze 





The continuity equation in the form 


Ou Ov OW 2.7) 
eo a 0 
Ox Jy Oz 


shows for a Cartesian coordinate system, the quantity of water 
entering an element must equal that leaving. The terms u, v 
and w are discharge velocity components of directions 


ay ard 


mammuespectfully. 


According to Darcy’s Law, expressed by equation (2.3), the 


component of discharge velocity can be expressed 





ah oh oh 
Lt = =k ae Y= -k —— , w= —| a) 
3 ay d= 
Sesser tuting into equation (2.7) 
—~k (ON/Ox) —kidh/dy) erage ol ? (ar) 
———— + 93——— + 3——— =-0 


Ox ov Oz 
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If k is constant for the element and then expressed ina 
typical two-dimensional form, equation (2.9) becomes the coamman 
Laplace equation. 


5 2 2 


The general conditon of flow analysis using flow nets is 
represented by Figttrre 2.11. Rules and good examples are 


presented in many textbooks (Cederogren, 1977s Harr, 1967). 


Free water surface 
(uppermost line of seepage) 





Flow 
lines 


Equipotential, 7 + 


lines 


Figure 2.11 General condition for flow lines and 
equipotential lines in flow net construction. 
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Dupuit Assumptions 
meee OS 


For unconfined seepage problems the engineer must calculate 
the location of the free surface. Many problems of flow have 
been approximated based on the Dupuit theory of unconfined 
flow. Dupuit (1893) assumed: (1) for small inclinations of 
the line of seepage, the streamlines can be taken as horizontal 
and equal potential lines vertical, and (2) the hydraulic 
gradient was equal to the slope of the free surface and 


imvarient with depth (Harr, 19642). 


Harr (1962) gives many example problems and has derived the 


2 ee 


condition of free surface elevations as shown by Figure 2.12. 


Figure 2.12 Simple free surface notation. 
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~ 5, 


In the case of Figure 2.12 the Dupuit formula 165 expressed as 


qg=k er) 


Many other reports and textbooks have presented charts and 
nomagraphs based on these assumptions. Later in this study a 
nhomoagraph for a simple trench drain with blanket drain 

SsOlution developed by the Federal Highway Administration (1980) 


Will be presented. 


Numerical Analysis 


Flow through porous media has been approximated by many authors 
USing ordinary non-linear differential equations and partial 
differential equations. Coupled with the speed of the computor, 
numerical solutions to many seepage problems have been 
developed. It should be pointed out that complex mathematical 
problems solved by numerical methods, are only as accurate as 
the math model can represent the conceptions and phenomina af 
seepage flow. In most cases, Simplification is the rule, and 
numerical models can only be treated as approximations to the 
actual occurances. Nevertheless, the number crunching 
capibilities and time savings involved gives the engineer much 


needed time for the true analysis of the seepage problem. 


The analytical and graphical methods mentioned previously are 
sometimes difficult and limited to flow systems which are 


ZT 





relativly simple. Im actual field conditions, the soils may be 
heterogeneous and anisotropic and boundary conditions very 
complex. These conditions make closed form methods very 
difficult and many simplifications would have to be made. 
However, with the numerical models, these complexities can be 
easily handled and many more check cases can be made in mtuch 


less time. 


Finite Difference Methods 


Many numerical computer models have been developed using finite 
difference methodology. Since the classic work of Feaceman and 
Rachford (1955), the ail industry has sponsored many studies on 
model development. A good historical background of the finite 
difference methods developed in seepage analysis is giving by 
Huyakorn (1983). Worth mentioning are the works of Cooley 
(1971) and Freeze (1971) in the development of other finite 
difference models. These models solve for the free surface and 


treat the saturated-unsaturated flow as an immobile air phase. 


Fimite Element 


The finite element method for solving subsurface flow has been 
widely used and researched be soils engineers and groundwater 
Hydrologists. The program used in this research and discussed 


in more detail in Chapter 7 is a finite element model. 
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Most finite element models are built on the Galerkin method. 
Some earlier finite element models for steady state and 
saturated flow were developed by Desai and Abel (1972), Conner 
and Brebbia (1977) and Zienkiewicz (1977). Finder and Frind 
(1972) used isoparametric quadrilateral elements with the 
Galerkin’s technique. Neuman (1973, 1974) used a four node 
quadrilateral element. A good reference for other works in the 
finite element medels that have Geen developed is reported in 


Huyakorn and Finder (1982) and Heijde (1985). 


Some significant works were those of Taylor and Hrown (1967) 
where a finite element model was developed for seepage through 
a dam solving for the free surface. The flow was only 
considerd for the saturated zone in this model. Freeze (1971) 
demonstrated that the conditions of the unsaturated zone may 
strongly influence the position of the free surface and results 
obtained not considering the capillary action for this case 
could be in error. Fredlund and Morgenstern (1976) developed 
improved constitutive equations for an unsaturated soil. Lam 
and Fredlund, etal (1984) presented additional theories and 
numerical development for further modeling techniques to 


consider the unsaturated zone in seepage analysis. 
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CHAFTER 3 


DIN)  eeeeteNT STUDY 


ael Finite Element Frogram Theory 


Introduction 


Finite element analysis of seepage problems can be described as 
follows (Huyakon and Finder, 1982): 

(1) A discrete number of interconnected nodes to define a 
series of sub-regions or elements. These can take the form of 
triangles, quadrilateral or rectangular shapes which are so 
spaced to define the geometric boundaries and soil types. 

(2) A Matrix expression is developed to relate flow 
Variables to each node of the elements. The mathematical 
expression may be of the form of either a variational or 
weighted function. A weighted residual function is used in the 
Galerkin method (Conner and Brebbia, 1977). 

(2) A set of algebraic equations that describe the entire 
global system 1s applied to assemble the element matrices. 
These equations are derived to consider compatible conditions 
at each node shared be each adjoining element. 

(4) Boundary conditions of flow are incorporated to reflect 
appropriate limits to the global matrix equations. 

(3) The resulting set of simultaneous algebraic equations 


are solved by some alogrithm method. Many algorithyms have 
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been used to efficiently take advantage of the banded and 
symmetric features of the coefficient matrix (Bathe and Wilson, 


ro7 Oo). 


The finite element method has great potential in seepage 
analysis problems, in that the complex geometries of flow 
regions can be easily constructed via nodes and elements. THis 
chapter will present the general theories and formulation of 


equations for the model used in this study. 


Frogram Theory 


A finite element program developed by Dr. Mustusfa Aral of the 
School of Civil Engineering, Georgia Institute of Technology, 

was used to study the flow and variables involved for 

subsurface drains as earlier described. This program was used 

to see effects of different parameters in the analysis performed. 
The program was operated on the Georgia Institute of Technology 
main-frame CDC CYBER computer. A complete user’s guide for the 
program 1s included as Appendix A-1l. ALSO example input and output 
used in this study is included by Appendix H-1. The formulation 
for the program presented here was summarized from a previous study 


(aieat. ocurm and Fultord, i981). 


This model only considers the flow condition of the free 
surface being the point where the water level is at atmospheric 


pressure. In other words, this model does not treat the 


oi 





unsaturated zone of the soil. For this study the flow continum 
is modeled to determine the quantity of flow at each drain and 
to determine the free surface location between each of the 
drains. This model also computed the velocities at each 


element for a good check. 


Equations of Flow 


Darcy’s Law is used as the basic relationship governing 
groundwater flow. Flow is assumed to be laminar and the 
inertial forces are negligible compared to viscous forces. 


Rewriting Darcy’s Law in the cartesian coordinate system: 


Vv =-K Va (3.1) 


where v is seepage velocity, Hh is piezometric head and FEF the 
hydraulic conductivity. In the two dimensional case, equation 


(3.1) can be written as 


dh 
Vv =- Kes — , iF 1,2 (S32) 


ox 


where Neg is a conductivity tensor. 
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The continuity equation for the flow in an element of soil can 
be expressed 


3 
Gen?) See Oo), t=zeii,e (35.3) 
1L ° 
Ox ot 
where 1 15 porosity of the soil and 0 18s the density of the 
fluid. Assuming the fluid and soil particles are 
incompressible, equation (3.2) and (3.3) can be used for the 


following expression in the two-dimensional form: 


= OO, i = 1,2 ae) 


Figure 2.1 can be refered to for notation of the key variables 


used 15 this study. 
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2 a Se a 


ds 


UOTIBION 9TQeTARPA YITM WaqIsSkS uTeApAspun TeoTdKky] 


Y 7 i" * doa 


ee ee ee - o— « 


———— —= — ee ee —Grem mute ss —-~--= 


ds ds ds 


l © Ofte my 


4H 
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Boundary COmnuL cLons 


For free surface problems, which is the case for this study, 
the exact location of the phreatic surface is not known at the 
beginning of the solution. Location of the free surface 15 an 
iterative process using two independent boundary conditions 
avallable at the free surface. The line of seepage at the free 
surface 1s the upper streamline of the flow domain. At every 
paint along this streamline the pore pressure is at atmospheric 
pressure and thereforethe piezometric head equals the elevation 
head. Alsa the flow across this boundary is zero because of it 


being a streamline. 


Considering the seepage domain as depicted by Figure 3.1, the 
steady state boundary conditions can be described as follows: 

(a) Impervious Boundary (AB): The program considers this 
boundary as a streamline where velocity to this boundary equals 
zero. 

(6b) Constant Head Boundaries (AC) and (EF): Comstant 
hydrostatic pressure is exerted on these boundaries where total 
piezometric head is constant, 1.8. equal to the water level in 
A plezometric tube. 

(Cc) FPReeatic Surface (CE), (FG), (HI), (JK) and (LM): This 
boundary 1s a combination of both cases (a) and (b) above. By 
iteration, described later, both conditions must be satisfied 


independently. 
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(d) Seepage Faces (EF), (GH), (IJ), (KL) and (MN): Seepage 
forces occur at these boundary faces where water exits through 
the soil continum, 1.@e. trench drains. Therfore streamlines of 


flow must be normal to these botundaries. 
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see Finite Element Formulation 


General Equations 


The first step in solving the seepage problem by the finite 
element method is to discretize the flow region into a finite 
number of subregions called elements. Using finite elements, 
equilibrium equations must be made to handle the geometries 
involved. This results in a set of Simultaneous algebraic 
equations rather than differential equations. Figure 3.2 
depicts the quadrilateral elements used for this program and 
shows how four six nodal triangles are formed internal to the 
program. A more detailed description is of finite element 


equations and formulations can be found in Zienkiewicz (1971), 


Desai (1972) and Aral (1974, 1976). 


ELEMENT CORNER 
NODES 


MID NODES 


Figure 3.2 Typical element and nodes 
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A finite element approximation to equation (3.5) can be written 
through a variational approach. Assuming Dirichlet and Neuman 
boundary conditions, and considering the two dimensional case 


af steady flow, the following equations can be expressed: 
h = oh on SS, a) 
and V = Sete oO ta Sn a OU on 35 7) 


where h is the specified Direchlet boundary condition in terms 
of piezometric head; Vv, is the velocity in the direction 
normal to the boundary; %, Sn peeves SCG iseGki on COSines of 
the normal with respect to x, and % § and > and 35 indicate 


respective boundaries of the domain. From equations (3.5) and 


(3.7), the following variational statement can be written: 


i (2 | i» 3 an\) 
ore mah —— + Sek. = hd, dx 
alr) Ze l Zz 
/ A Oh l 0 a IX. J a 


= Sy, ~ vs ShdsS (2.8) 
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Integrating by parts one obtains 


oh Oo Sh dh Oo dAh 
bs: —— ———— ClO Cl 
11 oy Ow ’ 0 x 0 x 


A 1 2 


| Ve 6hdS = 0 (3.9) 
55 


Since from earlier discussed boundary conditions as for this 


eeeem Studied, Vo =U. Thus, equation (3.9) can be further 


simplified by eliminating the boundary intreograls. 


Assuming equation (3.9) 1s written for a single element, matrix 
equations can be generated using an interpolation function for 
approximating Hh aver an element. The interpolation function 


used in this model is given as 


where h,*s are the unknown nodal values of the dependent 
ees ame, 165 the interpolating polynomial used to 


approximate h in an element. 
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For a typical six-nodal triangle element these polynomials are 


given as 
ay ee (een 1 ) 1 — tleea eae 8) 
and N. = 4E.E, eee Peed Gex1 >) 


where rr) » and b 4 are area cordinates defined by the 


following relations 


A 4 Cee) 
—— &, =n b,- — 
A 


where Al, A2 and Aj3z are the areas of the three subtriangles 
subtended by the point (F) and the corners, the index on (A) 


designating the opposite corner numbers, and Ais the area of 


the triangle (123) shown in Figure 2.3 


s(0,0,1) 





Figure 3.32 Area coordinates. 
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To simplify the resultant equations, the origin of the local 
axis which defines the local coordinates of the element, nodes 
are placed at the centroid of the element and the principle 
axis are inclined angle @ in the direction of local 


anisotropy. 


Substituting equation (3.19) into equation (2.9) the final form 
of the approximating equations in an element can be expressed. 
The stationary condition of equation (3.9) is obtained by 
equating the first partial derivatives of the function, with 


respect to the undetermined parameters, 3 TO mme mos 


Phys, = e = a ys a ON h. dx. ch: 
‘ 4 5 al 5 ; 5 - 299 - 7 fl es » 4 
ae a 2 
‘Ae 
Psd = 1,2,5;-.-4 (2.14) 


. : . = 
This results in a Six by six element stiffness matrix, S°, and 


e 
a six By one element load vector F. 
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The six by six element matrix equations of equation (3.14) is 


written 
aes) faa = eet GES 9. (3,15) 
where 
aN, ON, JN. ON 
[Ss] = eS dx. dx .. 
1 22 : ee 
Ke 1 1 2 2 
| = SRS s, PFS rae iG) 


The same procedure is taken for all elements. Then the total 
stiffness matrix is arrived at by combining the element 
stiffness matrices using structural assembly techniques 


CamemeWitcs, L971) 


For a grid of elements, a system Of banded simultaneous 


equations results: 


Bo 


ox 
ey 


, = Stes cane ) 


where ES] is the global matrix of cofficients which 

incorporates the properties of the materials and the geometry 3 
the elements; ths; is the vector of the unknown h*s at the 

nodes; and (Fi is the load vector. After the introduction of 
appropriate boundary conditions, this system of equations can 
now be solved By computations for the unknowns, hH, yielding the 


piezometric head distribution in the solution domain. 
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Free Surface Computation 


The flow domain must initially be established assuming an 
arbitrary free surtace, 1.@. the user initially estimates the 
free surface loacation. Al] boundary conditions are solved for 
the preset conditions of the problem. The nodes which are not 
fixed are allowed to move, where during the adjustment the 
interior nodes are also moved in order not to have distortian 
af the elements. Then additional iterations are made until all 
boundary conditions and user prescribed error correlations have 


been met for the configuration. 
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o.o Variables and Dimensional Analysis 


The variables investigated in this study are shown in Figure 
oe2l. For the analysis of seepage using the finite element 
program, variables were changed for each computer run. From 
each run data was generated to compare the changes occuring for 


each variable change. 


Only the two-dimensional flow condition was considered in this 
Study. For this flow condition, the building drainage system 
is i1n one cross-section of the building where each drain 
cOllects the seepage and collector pipes or drains take the 
drainage Dy gravity to a sump Or some exit drain. This is 


depicted by Figure 3.4. 
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co oe SUMP 


TRENCH DRAINS 





Figure 3.4 Underdrain System For Seepage Removal. 





The total discharge of the system, Ges can De express as the 
following function: 


= fe Dn aac. Saw 5. Sag 6 > cenit) 


a p? 
where Gy, = total discharge of all drains, i.e. the sum of 
discharges through each seepage face for the 
model, ft? /day/ft 
Hr = depth of saturated soil from the drain 


discharge elevation to the impervious 


boundary, ft 


H = fixed head of above the drains , ft 
dy = drain width, ft 
Sp = drain spacing, ft 


Ke = hydraulic conductivity of the soil, ft/day 

SFS = difference in elevation of the free surface to 
the water level of each drain, ft 

6 = slope of the impervious boundary with 


horizontal, ft/ft 


Dimensional analysis of equation (3.18) yields: 


Ceeveh 7) 


Since permeability is moved to the left side of equation (3.19) 
these variable is made proportional to the discharge. The 


Glevation of the free surface, Spos is an output variable, so 
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it 165 not given in equation (3.19. Also since the slope of the 
impervious boundary, 9, 15 already dimensionless it can be 


given as a Single term in FT/FT. 


The effects of these functional variables and other results are 
given in the form of dimensionless plots in Chapter 4. A 
discussion of the results 15 given in that chapter. Since 
considerable interest 15 in quantities of discharge and the 
location of the free surface, this analysis can be used to see 


trends and effects that each variable has on the seepage 


domain. 
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244 Frogram Verification 


To verify that the computer program used in this study is 
giving valid results, two other methods of solutions were used 
for a selected check case. Another micro computer finite 
element program, SEEF, and an appropriate analytical method 
based on Dupuit"*’s assumption, was used for the simple 
configuration of a drainage Blanket with a single pipe type 


“rT 
= 


drain (Duncan and Wong, 19385). This 18s shown in Figttre 2.35. 


afl 





Figure 3.5 Check case example of a blanket drain with single 
pipe trench drain. 

A grid was generated for the SEEF program to match the 
configuration and boundary conditions as by Figure 3.4. The 
total flow for this check case 15 calculated By the SEEF 
program as based on the Dupuit assumption of a fixed domain and 
fixed amount of head. Seep also tabulated the flow quantities 


for all nodal points with fixed heads. For this model, the 
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free surface converged on the free surface location with an error 


criteria of 9.16 feet. 


The location of the free surface as computed by the Georgia 
Tech model converged to the free surface Raving an error factor 
of only 90.001 feet, a much more accurate and sensitive 
convergence. However, the two models compared very close as 
far as the location of the free surface. Table 3.1 lists the 


elevations of the free surface nodes for the two models. 
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ieee >. 1 Elevations of the Free Surface Nodes 
Seer: VE I Ga Tech 
Model Model 
x coordinate y coordinate 
Q 73. OO Vee OU 

At) Vie a 7 eh 

70 67.07 C77 

BO 54.27 Sail 7 

Bo 62,03 63.84 

B75 60.39 sees) 7 

Ae. Se ot) 28,.00~ drain location 

pee 28. OO See) 

es 60,08 esol ha 

Oe wd Gia iG 62ers 

LOO Bun G2. 2 

106 63.00 San 2o 

Paes 51276 BB 

PAO 64.26 64.42 

150 64,42 564.62 
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To verify that flow quantities were accurate for the Ga Tech 
model, results were also compared with the results obtained 
from a Dupuit assumption solution, developed by the Federal 
Highway Administration’s, Highway Subdrainage Design Manual, 
(1980). Figure 3.6 was used to calculate the quantity of flow 


to the drain. The following calculations based on Figure =. 6 


are presented: 








b 1 Ls o 
— = = 0.017  , =) eae eo 
H, Sie ils we 
Pit — Ho? aes aa" ) 
is oe ee Time Ligure 3. o) 
2 42 


Ce) 





Sc@era! 2 N50 = ———— = 1.84 CFD 
ee ker 
q 1.84 Gq (H —- H )2 (30) 
= eT a ly a 
i O.2 K 2 (eee ye) EG i—-1) 
Qd 
=i) | lees PF eee) 6S UD aa 
E. 
therfore cs O.2011.44) = 2.29 CFD 


The results of this computation is compared with the results 


SSA 
Pratl 


of the two finite elements in Table 3. 
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Figure 3.6 Chart for Determining Flow Rate in Symmetrical 
Underdrains (after FHA, 1980) 
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Table 3.2 Discharge to the drain check cases, CFD 


Discharge, CFD 


SEEFsVFI Model eles 
Ga Tech Model ee 
FHA Nomograph Zn a? 


In summary, the Ga Tech model compares very well with other 
model and with the FHA nomograph for flow quantities. The Ga 
Tech model was found to more error criteria sensitive. The 
free surface location check also proved to compare very well 


with the SEEF model. 
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CHAFTER 4 
RESULTS AND DISCUSSION 
4.1 Introduction 


This chapter discusses the analysis and results of the two 
dimensional, steady state, free surface finite element study of 
Seepage to building underdrain systems. Details presented 
include model set-up such as boudary conditions, grid 
development and basic first assumptions. Output variables 
considered are total discharge and discharge at each drain, 


plus the free surface elevation between the drains. 


In stutdying sensitivity of the model for the input variables, 

it was found that the distance from the Dirichlet boundary 

face, or constant head boundary, Li, to the first drain was 
critical. Li was found to be critical in terms of output 
cCharcteristics, i.e. discharge and location of the free surface 
between drains. Therfore, two sets of results are presented for 


two Li distances of 400 feet and 80 feet. 
4.2 Model Set-up 


Figure 4.1 shows the basic grid that was developed for this 


Seepage problem. Using 150 elements and S21 nodes, the 


geometry of the mesh was Bee eg on a five drain system. The 
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Sides and bottoms of the drains were modeled as seepage faces. 
Using the two-dimensional model, the fifth drain is considered 
to be at the building centerline. Assuming symmetry, only 
one-half of the Building was analysed. The left most 

boundary, as shown in Figure 4.1, is a Dirichlet boundary 
condition which has the constant head required to allow for the 


seepage to enter the drains. 


Listed below are some basic assumptions for the study made for 
this seepage problem: 

(1) Steady state flow conditions 

(2) Two-dimensional analysis where the right most 
boundary is the building centerline and by symmetry the flows 
are only one-half of the total for the total width of the 
building. 

(3) Free surface nodes move to the computed phreatic 
level, 1.6. piezometric head equals atmospheric pressure at the 
upper most nodes 

(4) Fartial saturated flow is not considered 

(3) Evaporation 15 not considered 

(6) Retardation is not considered 

(7) Infiltration is not considered 

(8) The soil is homogeous being either isotropic or 


anisotropic. 
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4.2 Variables Used 


For this study, key variables were used in the finite element 
model to look at their effects on the flow of seepage for the 
building underdrain system. Two major design parameters 

necessary to model any seepage problem is the total discharge 
of the drain system and the location of the free surface. For 


this study, the following output variables are used: 


fee oreciarge Gr Seepage quantity ta each drain 
(cubic feet per day per linear foot) 
G7 = total discharge of all drains (same units) 


5 Fe = height of free surface above drain elevation. 


Discharge quantities are in units of flow quantity in time per 
unit length of the total flow continum. Therfore, flow for a 
total system should be multiplied be the length of the 


structure, 1.@. in the third dimension. 


Input variables used in this problem to study their effects oan 


seepage are listed below: 


Li = distance from the left most constant head 
boundary to the first drain (feet) 
E = hydraulic conductivity (feet/day) 


Evy = hydraulic conductivity in the vertical 


5/ 





direction (feet/day) 

Bh = hydraulic conductivity in the horizontal 
direction (feet/day) 

h = hydraulic head of flow continuum above the 
drain @levation (feet) 

He = depth of aquifer between the elevation or the 
drain and the impervious boundary (feet) 

Sp = spacing between drains (feet) 

Dw = drain width (feet) 


6 = slope of the impervious boundary (CA) 


Sweeter fFeference to the notation of the variables refer to 


sete. 


Figure 3.1 in Chapter :. 


For this finite element problem the units at feet and days wers 
usecd to give results for the output variables in cubic feet per 


day and feet per day for velocity. 


Hydraulic conductivity or permeability values used for this 
study were 09,0008, 0,02, 0,2, and 4 feet per day. mA osorl?s 
hydraulic conductivity having a value of G.,G002 feet per day 16 
equivalent to a fine grained sail such as a clay. Hydraulic 
conductivities are tabulated in Table 4.1 to Show cross 
reference ta soil types and permeabilities in units of 


centimeters per second. 





Table 4.1 Hydraulic Conductivity Values For Various Soil Types 


. cs | ; _-/ 
Fine grained soil O.9002 feet/day X10 cm/sec 
ae tae meetin: ’ er / 
rine grained soil 0,002 feet/day x#1O cm/sec 
es acer 
silt O.,O2 feet/dav 7X10 cm/sec 
sand, silt O,2 fFeet/day 7R1O > Sem sec 
_a/ 
clean sand 4 feet/day L.4x10 cm/sec 


The other input variables used were typically realistic values 
to meet typical field conditions. Listed below are some of 
this values, where for each computor run specific variables 
were changed and the others held constant. The intent was to 


see what effects that each variable change made on the model. 


Lio = 80°, 4007 and 1900° 
Hh = 10°, 20° and 40° 

igs = 20 reegomeand 15 
Soe la eos Aner aS" 

Dw = 1°, 3° and 4° 

oS) = OF nS... Wieeeande2o” 
RV =e COOS, O. oe 4 aad 15 
EA = O.2 


had =+ ad 


be ==) (OCG) & Q O,O2 3 4 
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4.4 Model Sensitivity 


The model was developed assuming a homogeneous isatropic media, 
except for a few cases where anisotropic conditions were studcdiec. 
The nodes of the top three rows of elements were allowed to move 
for iteration of the free surface. In all runs, the maximum 
nunber of iterations was Six and on an average three to four were 
found to be adequate. The model was mace to iterate to an 


accuracy of 0.001 feet during the free surface moves. 


In calibration of free surface models careful attention must be 
made to the moving nodes. Froblems exist when the elements becone 
less than zero in area or when the elements actual may become so 
distorted that they may become inverted. This was not the case 


for the mesh used in this study. 


In studying the flow quantities and velocities for the basic 
geametries and variable changes, all results seemed reasonable and 
observed trends appeared physically possible. However large 
variation of output was noticed by the change in distance from the 
constant head boundary to the first drain, Li. The first set of 
results is presented based on the influence length, Li, being 10 
times the maximum head difference, 1.8. 19 x 40° = 400°. Depth of 
aquifer for this case was held constant at 100 feet. Use of 10 
times the head differential is based on Bear’s (1979) 


recommendation. 
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Figure 4.2 depicts the moved free surface corner nodes for an 
influence length of 400 feet. For comparison, Figure 4.3 shows the 
moved free surface corner nodes for Li equal to B80 feet and a 


thickness of aquifer, Hr, @equal to S8 feet. 
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Figure 4.4 and 4.9 are presented ta show the dramatic influence 
that the distance from the constant head boundary from the first 
drain, Li, has on total discharge and the maximum free surface 
elevation, respectively. For further reference, the total 
discharge, q .i1s the sum of the drain discharges and the maximum 
free surface 1s the highest elevation of water above the drain 
@levation. Lengths between 200* and 400° appear to result in less 
change in flow quantity and free surface changes. This equates to 
the distance where the constant head houndary away from the drains 
begins to cause the least effect on discharge and free surface 
changes 15 approximately 15 to 20 times the head differential. 
This perhaps is good evidence that the actual field conditons for 
each particular application of groundwater movement should 


carefully be looked at in terms of this distance. 
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Figure 4.4 Distance from the constant head boundary to the tirst 


drain vs. 


h = 20°, 
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4.5 Results 


This section will present the computer results in the form of 
dimensionless plots and other nomographs to show the effects 
that each variable change has on the seepage problen. Two 

sets of plots are given where the first set is for the constant 
head boundary distance, Li, equals 400 feet and the depth of 
rock, Hr, equals 100 feet. The second set of plots are 
essentially the same as the first except Li equals 80 feet and 


the Hr equals 3S feet. For variable notation refer to Figure 


Ls 
Drain Discharge (Li = 400° and Hr = 100°") 


Figure 4.6 shows the influence of drain width, Dw, for head 
differentials, h, of 10, 20 and 40 feet on total discharge of 
the drains. This shows that drain width is not particularly 
critical for low head differences. Tt 15 reasonable to assume 
that total discharge could approach a constant for drain widths 
greater that four feet and head less than 2O feet. Although 
total discharge increases with drain width, this plot does not 
show what drain width is the most effective or most efficient. 


However, the influence of head can be easily observed 
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Figure 4.6 Total drain discharge vs. drain width with head 


differentials. 
Pee. Hr = 100we Sp = Soy = O.e ft/day 
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Figure 4.7 depicts the results of varying anisotropic soils 
With varying heads. In this analysis, vertical hydraulic 
conductivity varies while horizontal conductivity was held 
constant. Figure 4.7 shows that with increase in vertica 


hydraulic conductivity discharge increases. This appears to 


satisfy theory , that with less resistance flow increases. 
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Figure 4.7 Total drain discharge vs. anisotropy with varying 
head differential. 
Li = 400°, Hr = 100°, Sp = 15°, Dw = 2’, Fh = Of ft/day 





Figures 4.8 through 4.12 are graphs which show the discharge 


Quantities for each of the five drains. 
effects of depth of impervious layer, 
discharge for each of the drains. 


increased depth discharge increases. 
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Figure 4.9 shows the effects of head differential on the 
discharge of each of the drains. It is interesting to see that 
though the first drain is the most critical drain, i. @. a 

large amount of the total flow is collected in this drain, the 
remaining drains collect discharge decreasing linearly from 

the left. This is possiblly due to the influence of the 
constant head boundary distance, Li, being 400 feet away. This 
Shows that for the drains, flow quantity increases linearly with 


increase in head. 
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Figure 4.9 Drain discharge vs. head di erential. 
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Figure 4.10 depicts the results of Comparing the drain 
discharge versus the spacing between drains. It is of 
interest to see that drain number one for spacing equal to 25 
feet, collects more seepage than the other discharges for the 
other spacings. However, at drains 2 through S&S this changes 
and the discharge for 15 feet spacing is greatest. This will 


be analyzed later in the dimensional analysis. 
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Figure 4.10 Drain discharge vs. drain spacing. 
bi = 4007, He = 1007, h = 20°, Dw = 3°, EK = 0.2 £t/day 
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Figure 4.11 shows the effects of anisotropy of the soil 


matrix on drain discharge. It is of interest to see that 
the discharge for Ev/EhW equal one 1s much greater than the 


other drain discharges for Ev/EhH less than one. 
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Figure 4.11 Drain discharge vs. anisotropy. 
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Figure 4.12 depicts data generated for the condition of a 
sloping impervious layer. The slope of the boundary is given 
as a percentage of Ax/Ay. The depth of the layer is held 
canstant at drain one and the layer slopes clockwise from the 


positive x horizontal axis. 


Though actual quantities of discharge are only changing in very 
small ammounts, the trend is seen that for slopes greater than 
wm the discharge increases in drains 3 through 3. However at 
some slope between SZ and 10% the discharge decreases. It is 


also seen that the discharge of drain five is the least effected 


by any slope change, in terms of discharge quantity. 
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Free Surface Location (Li = 400° and Hr = 100") 


Figure 4.12 shows the influence of the drain width on the 
maximum free surface height, sero In all cases this occured 
between drains one and two. This plot shows how the 


relationship is linear and that the height of the free surface 


1S not greatly affected by the drain width. 
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Figure 4.14 depicts the influence of drain spacing and head oan 


These results show that 


the maximum free surface height. 


freer 


increased drain spacing causes an increase in the 


cor 


Therefore one should consider the smaller spacing 


surface. 
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Figure 4.14 Maximum free surface height vs. 


with varying head. 
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Figure 4.15 shows the relationship of maximum free surface 
height, -FS ax! as a fuction of head, h, and anisotropy, Kv/Eh. 
The same trend is seen here as noticed with discharge for 


anisotropy. that with increased vertical hydraulic 


canductivity the height of the free surface increases. 
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Figure 4.15 Maximum free surface height vs. anisotropy with 


varying head. 
Li = 400, Hr = 100°, Sp = 25°, Dw = 2° 
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In Figures 4.16 through 4.20, the free surface corner nodes are 
plotted for the different cases of analysis. To see the 
effects of the free surface changes versus the variable changes 
only nodes between the drains are plotted. In the sence of 
design, the maximum height of the free surface would be of the 
most concern, which 15 in between the first two drains. 
However, it 1s of interest to see the total free surface from a 


standpoint of total cross section of the building. 


In Figure 4.16 the effect of thickness of aquifer on location 

of the free surface is shown. It is of interest to see, as in 

the case for discharge shown by Figure 4.8, that for shallow depths 
Qf impervious layer, Hr = 20°, that very little rise in the 


free surface occurs in between drains four and five. 


Figure 4.17 shows the effects of head differential on the 
location of the free surface nodes. It 16 of interest to see 
for cases of head equal to 20 and 10 feet the free surface 
nodes are approaching gentle mounds of water in between drains. 
However for the head equal to 40 feet the free surface peaks 
rather abruptly between the first two drains but resumes the 


gentle fall in free surface between the other drains. 


Figure 4.18 shows the effects of drain spacing on the location 
Of the free surface nodes. This shows as spacing is increased 


the free surface increases. 
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Peorgere +.27, tmeserrects Of the sloping aquifer on location 
of the free surface nodes is shown. The same trend is seen 
here as in discharge, shown by Figure 4.12, that for slopes 
greater than 10% the free surface location drops. The change 
in free surface nodes is very small for the cases of slope 


equal to SA and 10OZ. 
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Figure 4.19 Free Surface corner nodes with drains with effects of 


slope of impervious boundary. 
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Dimensionsal Analysis (Li = 400° and Hr = 100°) 
Piette Se TT LS 


Figures 4.21 through 4.24 are plots of the data generated Dy a 
dimensional analysis as discussed in Chapter 3. The discharge 
values given are unit discharge quantities and each case is for 
the particular variables noted in each figure. Unit discharge 
is defined as the sum of the total drain discharges for each 

of the drains per unit width in the = axis or the third 


dimension. 


In Figure 4.21 the effect of drain width is analyzed. This 
plot shows that with increased drain width discharge increases, 
but with increased head the discharge may increase at a much 
faster rate. From this analysis it would appear that the 
greatest increases of flow rate would be for the drain widths 
between one and two feet wide. The following analysis 15 


performed: 


Drain Widths 





head 1 = as 4 
ane) QO, O24 iO cys 0.065 a7 S 
increase —- 0.019 Q.Q12 QO, O09 
20 Oi Og ees 0,070 0.986 
increase —Y 0.02 Ons O.O1L16 
40 0,036 Oa / O,O73 QO.09 
increase —e, 0.021 0,018 geal s 
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This analysis shows the increases of the flow fumetieonm for the 
different head increases as shown. The maximum increase 1s 
seen for the values between drain widths of one and two. 
Therfore the most effective increase in flow occurs for drain 
widths between 1 and 2 feet. Another relationship could be 


= 


stated, that for Dw/h ratios less than 0.2 the flow quantity 


increases at a rapid rate. 
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In Figure 4.22 the effects of thickness of aquifer is shown for 
the dimensional snalysis- This plot concludes that for Hr/h 
ratios greater than 2 or 2 the net increase in flow is 
approaching a linear relationship. This states that with 
increasing thickness of the aquifer the dicharge is expected to 
increase very near a linear increase and for very shallow 


thicknesses the flow becomes very small. 
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Figure 4.22 Dimensional Analysis, Depths of impervious layer. 
li = 400°, Sp = 25°, Dw = 32", K = O.2 ft/day 
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The effects of drain spacing is presented in Figure 4.23. The 
flow function for this analysis 15 treated differently for the 
spacing changes. As the spacing was increased for each of the 
computer runs the flow domain increased accordingly to 
accommodate the five drains at the appropriate spacing. 
Therfore to gain understanding of the effects of flow 
relationships with respect to drain spacing, each term was 
normalized by the corresponding spacing between the drains for 
each head increase. This analysis shows for d@creased spacing 


discharge increases and as head increases so the discharge. 
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Figure 4.25 Dimensional Analysis, Drain spacing. 

Li = 400, Hr = 100°, Dw = 2°, K = 0.2 ft/day 

Figure 4.24 shows the dimensional analysis for effects of the 
sloping aquifer. For this particular flow matrix as modeled, 
a definate break in flow increase is noticed. The peak 


discharge function is at slope, 9 equal 0.05. 
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Figure 4.24 Dimensional Analysis, Slope of aquifer. 
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Results For Li = 80° and Hr = 58’ 


The following plots and nomographs are the results for the 
computor runs where distance from the constant head boundary to 
the first drain, Li, equals 80 feet and the thickness of the 
aquifer equals 38 feet. The plots are essentially the same as 
for the previous case except these variables are held 


constant. 

Drain Discharge ( Li = 80* and Hr = 5g’ ) 

Figure 4.25 shows the effects of drain width on total drain 
discharge. The results are generally the same as shown in 


Figure 4.6 except discharge quantities range from = to 4 times 


greater. 
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Figure 4.25 Total drain discharge vs. drain width with varying 
head. 
mee eso. Hr = 58°, Sp = 25°, K = O12 Ft/day 
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The results of anisotropy of the soil continum is shown by 
Figure 4.26 for the short Li distance of 80 feet. The trend is 
slightly different than for the case for Li equal 400 feet as 
shown by Figure 4.7. For the largest head value of 40 feet the 
discharge decreases just slightly with decrease in vertical 
hydraulic conductivity, but for a Ev/Eh ratio less than 9.1 the 
discharge decrease is very dramatic. However for the low head 
Value of 10 feet the general decreasing trend is prevelant as 


the case for Figure 4.7. 
For this case the flow quantities have increased approximately 


4 to 3S times. This shows the large effects that the distance 


of the constant head boundary has on the seepage quantities. 
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Figures 4.27 and 4.28 depict the same general trend of maximum 
flow at drain one and a general decrease in flow tor each 
consectutive drain is noticed for effects of thickness of 
aquifier and head respectively. Flow quantities are again 2 to 
4 times greater for this case than for the earlier case as 
shown by Figures 4.8 and 4.9. For this case ot shorter 
distance to the constant head boundary equal BO feet, the 
mid-peak discharges at the third drain is not noticed as in 
Figure 4.9. In Figure 4.28 it is of interest to see that the 
largest head value equal 409 feet has such a large effect on the 
first drain collecting a large portion of the total PP eharce 


of the system. 
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The effects of drain spacing for distance of the constant head 
Poem to the first “drain, Li, equal SO feet, 1s shown by 
Figure 4.29. The same trends are seen for this case as that 
shown by Figure 4.10. Flow increase is somewhat less in this 
case compared to the earlier case of Li equal 4090 feet. Flow 


quantities appear to range from 2 to 2 times as great. 
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Figure 4.29 Drain discharge vs. drain spacing. 
li = 80°, Hr = S38’, bh = 2O°, Dw = 2°, KF = O12 Ft/day 
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In Figure 4.30 the effects of slope of the impervious boundary, 
G » is shown for drain discharge. A significant change is 
noticed for this plot when compared to the earlier plot in 
Figure 4.12 for sloping aquifer. For Li equal to 400 feet, the 
peak discharge occured for a slope between 3S and 10 percent. 

In this case, where Li equals 80 feet, the peak flow is for the 
case of the steepest slope equal to £9 percent. For drain one, 
the flow is 8 times greater than the flow for drain one as 
shown by Figure 4.12, Li = 400 feet. The flow is 3 times 


greater in drain five in this case than the flow in drain five 


for the case of Li = 400 feet. 
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ee eee eee ones Sh eee end Ar = 58’ 


Figures 4.21 through 4.26 show the results of the free surface 
location plots for the case of Li equal 80 feet and Hr equal S38 
feet. The same general trends are shown for this plots as 
shown by the earlier case of Li equal 400 feet and Hr equal 100 
feet. In general, the free surface has moved up approximately 
three times in height. An exception to the earlier case of is 
shown in Figure 4.26. This figure shows the location of the 
free surface nodes with the effects of slope of the impervious 
boundary. It is samewhat different in this case than in the 
earlier case as shown in Figure 4.19 for Li equal 400 feet. rn 
this case the free surface continues to increase in height as 


the slope of the impervious boundary increases. 
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Figure 4.21 
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Figure 4.32 Maximum free surface height vs. frain spacing 
with head varying. 
meee. He = 58°, Dw = 2’, K = 0.2 ft/day 
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Dimensional Analysis ( Li = 80° and Hr = SS" ) 


The effects of drain width as shown by the dimensional analysis 
is shown by Figure 4.27. The general trend of the discharge 
increasing with drain width is seen in the data. When compared 
to the earlier case for Li equal to 400 feet, it is seen that 
total discharge is greatly effected, by a factor of three 
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Figure 4.37 Dimensional Analysis, drain width. 
Peo. rh = 28°, Sp = 25°, KF = 0.2 ft/day 
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Figure 4.28 shows the dimensional analysis results for 
thichness of the aquifer. In this case, Li equal 80 feet, the 
relationship is more linear than in the case of Li equal 400 
feet as shown in Figure 4.22. This results in the conclusion 
that the effects of Hr has a pronounce influence on the flow 
quantities when coupled with the distance Li. For a shorter 
distance of Li, the effects of depths of Hr is more crictical 
and more of a linear relationship. However when the distance 
from the constant head boundary is larger, as in the earlier 
case, the effects of Hr are not as effective when Hr is very 


small. 
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Figure 4.39 shows the results of the dimensional analysis for 
drain spacing. Again the ordinate is normalized with respect 
to the area subtended by each flow continum for each case of 
spacing. Flow quantity would equate to an increase of 
approximately S to 6 times from this case of Li equal 980 feet 


an ae 
versus Li equal 400 as shown in Figure 4.<.. 


— 
> = 
5 = $s 
= 

— S a — OP RTTO TEE Wat ee ear mS FS eT Fl TL TT TT TS A A I yA Re 


= = - 2 ow wm 


begyeat 
3 oh 


i 


Y 


a? 


AX 
OT 


haet 
rent 


igts? 
) 


pity 





til 
inti 
"| 


Hh 


: : . 
: : . 
CT I A AR A PS Me St EE AD nS SY tL A Se aC See SS of hemi repr 
= : * : 


H0/;, 


*) 
Vas 
a 
yj 
t 
ae 


( 


? 
2 > 


- 


t "| 
ae | 
{ 


Figure 4.39 Dimensional Analysis, drain spacing. 
Li = 80°, Hr = 58°, Dw = 2°, FF Oe #t/day 
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Figure 4.490 shows the dimensional analysis results for the 

sloping impervious boundary. For this case, Pimedtial Steer, 
as depicted in the earlier plots for sloping aquifer for drain 
dischrge and free surface location, the trend is seen that for 


and increase in slope the discharge increases. 
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4.6 Summary 


This chapter has included the results of many computer runs for 
the analysis of a Building underdrain system. The importance 

of model sensitivity must be stressed, in that key variables may 
effect the results very significantly. In this model, it was 
found that the most crictical variables affecting the total 
discharges of the drains and location of the free surface was 
the distance from the constant head boundary to the drain 


system, Li, and thickness of the aquifer, Hr. 


The results obtained appear to be reasonable and within limits 
of possible field occurances. However, using these results for 
quantifying flows and the free surface locations requires a 
very careful study and Background investigation of all 
geometric configurations and soil conditions. It 15 intended 
that this study has presented the results in a manner that one 
can observe the effects of the variables studied on such a 


underdrain system. 


In summary, the following observations are discussed for the 


variables used in this analysis. 


Constant Head Boundary Distance, Li 


As distance from the constant head boundary to the system 
Lo 





increases the flow quantities increase for the drain system. 
In terms of a correlation with hydraulic head differential 
above the drain, as used by other authors, this model study 
shows that for a distance of approximately 15 times the head, 
discharge appears to decrease. In terms of gradient, this 
equates to a l over 15 ratio or 9.9067 gradient. A slightly 
higher distance was observed in the analysis of this influence 


to location of the free surface. 


Depth of Aquiter, Hr 


The results show that with increased aquifer thickness the 
drain discharge increases, increasing very near to a linear 
relationship. It was also observed that for the very shallow 
thicknesses, the last drains collected very little flow. The 
height of the free surface is affected in the same fashion. 
With increased thickness of aquifer the height of the free 


surface increases between drains. 


Drain Spacing, Sp 


Drain spacing effects on discharge was greatly affected by the 
distance of the system from the constant head boundary. A 
substantial increase of 3 to 6 times the discharge was observed 
for the two distances used in the results presented. It was 
observed that as head increases, discharge increases for all 


drain spacings; and as drain spacing increases, discharge 


Aah? 





decreases. It was also observed that from the dimensional 
analysis, the most efficient drain spacing 1s for a spacing 

over head ratio of less than 2. The location of the free surface 
was observed to increase between drains as drain spacing is 


increased. 


Drain Width, Dw 


With increased drain width the discharge increases tor all 
other variables constant. However, there must be a trade off 
for efficiency of the complete system. It appears that a drain 
width between 1 to 2 feet is the most effective width. This 1s 
also a plus from a standpoint that lateral seepage can be 
easily drained into thinner trenches or deeper constructed 


drains. 


Slope of Impervious Boundary, © 


For the distance of the constant head boundary to the drain 
system equal to 890 feet, the drain discharge was observed to 
increase with increased slope. However, for the distance, Li 
equal to 400 feet, the maximum discharge occurs at a slope 
between 2 to 10 percent. This occurance 15 not fully 

understood and perhaps should be further studied to develope 

some reasoning. This requires that careful consideration of 
specific boundary geometries be used in model studies for seepage 


analysis. 
uv TS 





CHAFTER 3S 


CONCLUSSION 


In this day and age, more and more sites are chosen for 
buildings and engineered structures that are less desirable to 
work with from a standpoint of site conditions, poor soils and 
water problems. A less desirable site may easily be where the 
groundwater table is very near to the surface. Therfore the 
building’s design must incorporate some thought into how the 
groundwater will be handled and what effects it may have on the 


structure. 


Engineers faced with the analysis of seepage for buildings to 
be built below the groundwater table must consider such 
conditions of uplift pressures, critical velocities, seepage 
control of the water itself or particle migration. Hopefully 
the design will enhance stability of the structure and prolong 


the usefulness and life of the structure. 


In performing the literature review on this topic, it was found 
that many authors and researchers are more and more relying on 
and developing computor models to perform the analysis of 
complex seepage problems. With the use of general seepage 


theory and mathematical representation of the seepage theories, 
ie 





the computer programs can quickly and efficiently aid the 

endgineer in studying the flow conditions. As in this study, 
the analysis of one particular seepage problem using a finite 
element program made it possible to study many key variables 


and develope many relationships. 


For this study, underdrain systems for buildings built below 
the groundwater table were researched. By using a finite 

element program, and by developing a mesh for this particular 
flow matrix, many variables of the analysis were studied. Tt 
was found that the model could easily be used to efficiently 
and accurately study these variables. Of major concern, the 


dicharge of the system and location of the free surface was 


studied in detail versus the various varible analysis. 


This research has attempted to present general seepage theory 
and useful criteria necessary to analyzing this particular 
seepage problem. It has also attempted to use state of the art 
techniques, 1.@e. finite element models, to analyze the 
Variables of the underdrain system design. In conclusion, the 
author must urge that in using such models that the basic 
assumptions of the model be thoroughly understood due to the 
somewhat simplifications of the some model operations, and that 
the engineer obtain as much information and data for the soil 


parameters and fluid boundary conditions which may exist. 
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AFFENDIX A-1i 


USER’? S MANUAL FOR 


ARAL SEEFAGE FROGRAM 


l. INTRODUCTION: 


This seepage program 1s a finite element model for solving 
unconfined or confined two dimensional or axisymetric seepage 
problems. For unconfined flow analysis the model iterates by 


moving the free surface nodes to the phreatic level. 


In general, the program allows for the porous soil mass to be 
divided into discrete elements and nodes using either a vertical 
plane model (two-dimensional,x-y) or axisymmetric grid structure. 
Then applying known boundary conditions, the program 
Simultaneously solves for each element, the continuity equation 


which satisfies Darcy’s Law. 


Figure 1 shows a very simple case of a two dimensional unconfined 
Steady state condition which can serve as an example. Pima ceinn 
example the free surface location is unknown and the user must 
begin the process by initiallizing the input by approximating the 
location. Then the program will iterate to the free surface level 
until the free surface converges to satisfy the equations of 


flow. 





Figure 1 


EP 


IMPERVIOUS BOUNDARY 


Example Grid for Unconfined Flow Continum. 





Il. FROGRAM EXECUTION AND ANALYSIS: 


The user has the option of performing a two dimensional problem or 
an axisymmetric problem. The two dimensional case simply shown by 
figure 1 will serve as an example to explain the parameters. An 
axisymetric application might be a well source problem looking at 
stratified soils having varying permeabilities at different 


elevations through the porous media. 
III. BOUNDARY CONDITIONS OF FLOW: 


The potential, Q, of flow must be defined as some directional 


function 


pi = f (isynz4t) 


A Dirichlet boundary condition is where the piezometric head, h iz 


constant, or 


h = f (Xsyszst) 


is solved for at a particular point. This boundary condition 
occurs whenever the flow domain is adjacent to an open body of 
fluid. Since the piezometric head is constant at the open body of 


water, the piezometric head of the porous media at the interface 
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Will be constant. This boundary condition denotes a non- zero 


flux or flow is to be considered possible. 


Another boundary condition for flow potential is the Neuman 
boundaries. This 185 the case where flow cannot move through the 
porous media at a Specific point defined, 1.€@. an impervious lay= 


of stratum. 


IV. FINITE ELEMENT GRID NOTATION: 


The porous soil mass to be analysed must be divided into a 


grid of elements and nodes. The program uses any shape second-order 


isoparametric element, shown By figure 2 as a Quadralateral 


element containing four corner nodes and four midpoint nodes. 


-CORNER NODES 


MID-NODES 





Figure 2 Element and Node Notation. 





The nodes must be numbered from some beginning corner, preferably 
the number one node should be at the bottom left corner of the 
mesh. Then node numbering increases to the end of each column, 
beginning at the bottom of each column throughout the mesh. One 
can see the numbering scheme by the examples presented. Elements 


are mumbered the same as the nodes. 


The program automatically generates nodes at equal spacing between 
corner nodes which are defined by the user. This feature 

requires that careful attention be made to the node numbering scheme 
once corner nodes have been established. This feature is good for 


creating finer mesh grids in particular areas of interest. 


Rase elements are used to create stationary reference elements 
typically at the bottom of the mesh or away from the free surface 
elements. These elements are also used to define changes in 


Hydraulic conductivity in the x and y directions(kx and Ey). 


The free surface or phreatic surface of a media is defined by the 
user as a first approximation of it*s actual location. In the case 
Of a cutoff wall, shown in Figure 3, there is no free surface and 
the seepage model solves for the flow characteristics of each 
element and node. However, in the case of a dam as shown in 

Figure 3, the free surface is unknown and the program will solve 


for the phreatic surface where the fluid pressure is at 
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atmospheric pressure. In confined aquifers the water pressures 


would possible be higher than atmospheric pressure. 


Seepage faces are defined by the user to denote a boundary or 
where seepage is .to be allowed. Here flow is qunatified and the 


model solves for the velocity and discharge at this face. 





CUTOFF WALL DAM CROSS-SECTION 


Figure 2 Example of confined flow and unconfined flaw. 


i 





Fermeablility of the porous media often depends on the direction 
Of flow. For the condition of soil matrix being anisotropic the 
velocities and gradients are generally not parallel. For this 


condition the velocities can be written in the Darcy Law form: 


h 
Vs = =e > 3 


Ox 


Oh 
Vy yy 
ay 


Oh 


oz 


Normally the hydraulic conductivity of teh media is similar in the 
y direction. Ther fore treatment of this normal condition requires 
use of only two direcetions. For this simplification, the program 
uses the notation of the two dimensional case. Figure 4 shows the 
notation that should be used for describing the hydraulic 


conductivity conditions. 


KY 
Xt 


ALP 


Figure 4 Hydraulic Conductivity Notation. 
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Via INFUT VARTABLES AND DATA INFUT 


Input can be given in unformated notation where each data input is 
seperated by a comma or space. Careful attention, however, must 
be made to use the apppropriate number input data per line as 
described by this guide. 


Title and Type Froblem Card. 


aoe 


DIMEN 


FERR 


The title can be up to 72 characters. 


Flow problem e problem is a two dimensional analysis 
input 2 and for an axisymetric problem input 1. 


Flane Flow Froblem = 


Axisymetric Froblem 1 


Having performed a sensitivity of the model or decided 
on an error critiria, input the error required. 


Control cards 


NNF'C 


NELEMC 


NF ROC 


NFFS 


Total number of corner nodal points. 

These nodes are used to define corners of elements 
which are defined By input coordinates. Nodes that 
are not listed will be generated between two corner 
nodes for each vertical column or row deliniation. 


Total number of base elements. 

Rase elements are typically at the bottom of a series 
Of elements that change the hydraulic conductivity of 
the elements above these base elements. These elements 
are established as reference elements for the program 
to describe permeability of each element. 


Total number of Neuman Boundary Faces. 

A Neuman face is defined as a line of nodes that may 
create a barrior or impervious wall for the flow. Lf 
two corner nodes define a line of several nodes that 
makes an impervious barrier, this is one Neuman 
Roundary face. 


Total number of corner nodes on the free surface. 

The user must specify particular corner nodes which 
reflect the first approximation of the unknown free 
surface. The program uses these nodes and the midpoint 
nodes to converge to the free surface. 
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I T Gav Number of iterations desired. 
This 1s a safety feature for possible input errors and 
computor time limitation. 


NF EDC Total number of Dirichlet Boundary faces. 
Flow is allowed at these faces. Do not count the mid 
point nodes as part of the face. 


NSFFAC Total number of seepage faces with Dirichlet 
boundaries. For each face defined by two corner no jes 
for the NFRDC input variable a NSFFAC value should be 
given. AS in 


NFRDC = 2 NSFFAC = 1 
NFRDC = 3 NSFFAC = 2 
[TIME Time variable. Input © = time independent. 


This is used for transient computation and is not 
used im steady state computations. 
Corner Node Location Cards 
Corner node number (n) 
x(nd «x CcoGardinate of node (n) 
y(n) y coordinate of node (n) 
Note: There should be NNFC number of these cards. Use 


appropriate dimensions, i.e. feet, meters which will be carried 
throughout the remaining input variables. 


NFMIS Enter a "1" or "©" where 
1 - additional nodes are to be generated after this 
node. 


QO —- no nodes will be generated after this node. 
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Hase Reference Elements Identification Cards 


NOD 


NMIS 


Xk CN) 


YEON) 


ALF (N) 


For each base element(N), give a counter-clockwise 
numbering tor each node in the element. 


5 i 20 
a IS 
| tie 18 


Total number of elements above the base element before 
the next base element. 


Hydraulic conductivity of the elements defined by the 
base element in the x direction. 


Hydraulic conductivity in the y direction. 


Angle of possible stratification or plane of soil 
matrix with the horizontal. 


y 
KS 
Kt 
ALP 


Note: There should be NELEMC number of these cards provided. 


Neuman Houndary Condition Cards 


NSTART 


a= 


RVAL 


NESAME 


The beginning or lowest corner node number of each 
Neuman Houndary defined. 


The Neuman Boundary value that defines the known 
potential height or water level measured from the 
datum reference line of the coordinate system. 


Total number of boundary faces having same Neuman 


boundary condition. Each face is a line made by two 
corner nodes of the same Neuman boundary. 


A-1l 





NINC The number of nodal number increases or numerical 
difference to the next corner node from node NSTART. 
Do not count mid-nodes only corner nodes. 


Notes There should be NFBOC cards provided for this set. 


Dirichlet Houndary Condition Cards 


NSTAR The beginning or lowest corner node number which 
defines each Dirichlet boundary. 


DEVAL The Dirichlet boundary value that defines the potential 
or water level height measured from the data reference 
elevation of the coordinate system. 


NDSAME Total number of boundary faces of this Dirichlet 
boundary condition. Each face is counted between 
corner nodes defining the boundary or number of 
elements. 


NDINC The number of modal number increases or numerical 
difference to the next corner node from node NSTAR. 


NMID The number of nodal number increases from the NSTAR to 
the next mid-node. 


Note: There should be NFBDC cards provided for this set. 


Free Surface Nodal Cards (Maximum 16 pairs per card) 
NFFSA Free surface corner node number. 


NF ROT Node number of the bottom of the movable section of 
the grid that defines the limits of moving modes for 
NFFSA. 


Note: There should be NFFS pairs of node numbers given for a 
maximum of 16 pairs per card. Use as many cards as required. The 
bottom node (NFBOT) defines the bottom of a seepage face 
corresponding to the grid geometry and corner node defining the 
free surface. The program uses this band of nodes and elements to 
converge on the free surface. 
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Mid-point Free Surface Nodal Cards. 
NFFSM The mid-point nodal numbers between corner nodes 
defining the free surface (NFFS —- 1). 
Seepage Face Cards 
NSFC Reginning node number of the defined seepage face. 


NSSAME Number of seepage faces defined by elements for this 
seepage face NSFC. 


NSINC Node numbering increase between nodes defined by 
NSSAME. 


Note: There should be NSFFAC cards provided for this set. 


Statement and End Card 


End of Froblem Statement should be less than 72 characters and 
end with a period. 


NDIMEN Input a "3" to end the problem. 


ERR Specify any error as for FERR. 
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UNCONFINED FLOW EXAMPLE FRORLEM 


a 


QYRICHLET Bou PARY Cak/TTIOKS 





| 
r 


BASE ROFECREICE 


ELEM ELST S 

Title and Type Froblem Card 

TITLE - SILTY SAND DAM 

DIMEN - 2 (two dimensional ) 

FERR - 0.01 (error criterion is 9.01) 

Control Cards 

NMFC - 12 (corner nodes are 18, S32 35, 45, 32, 62, 69, 79, 86 and 

Oa 

NELEMC- & ¢( five base elements —- shaded) 

NFBOC - © ( no Neuman Boundaries) 

NFFS - 6 ( & nodal points on the free surface - 11, 28, 45, 6&2, 

79 and 96) 

ITGIV -—- 3O ¢( maximum number of iterations) 

NFBDC ~— 2 (There are two faces with Dirichlet boundary conditions, 
the left side and the right side of the mesh which 
forms the dam) 

NSFFAC- 1 (There is one seepage face —- on the right side where 
flow is to be considered.) 

ITIME -—- © (Time independent for steady state conditions) 


A-14 





——— Ieee Ch l(c 


Corner Node Location Cards 
foe). 1 


Where, the first node is number 1 having an x coordinate of © and 
a y coordinate of QO. Generation of additional nodes will follow. 


fee 2 10,0 

Where, node 11 is at the end of the first column at 
coordinates O and 110 and no generation is required for other 
nodes in this column. 


oO, 0,1 
Bono, 79,0) 
ade 400,090, 1 
45,400,88,0 
eeu), O,1 
ee 0), 77,0 
Been, 0,1 
ree. 66, 0 
mem OO, 0,1 
Toren Oo), Jo, 0 


Rase Reference Elements Identification Cards 

meee) 7, 20,15, 5,2,4,0.002,0.002,9 

Counter clockwise numbering of nodes for each element and the 
permeability, 90.002 ft/day for this base element and 4 elements 
above the base element. 

Pee y OO, 37, 50, 20,19,4,9.002,9.002,9 

te boy ae Soe 45 475275 SO, 4,0.002,9,.002,0 
Gewor.O7./70,/71,64, 04, 795,4,0.002,0. 002,0 

eee, 7, 0,51, 71,70,4,9.002,9.002,9 


Neuman Houndary Condition Cards 


None specified, so no cards will be provided 
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Dirichlet Houndary Condition Cards 


moe, J, oy lL 

fees, J, 2. 1 

The first node of the first Dirichlet Boundary is node 1 and the 
DEVAL value is 1190 feet. There are S elements on this face that 
makes this Dirichlet boundary condition, NDSAME = 3S. NDINC equals 
2 for the nodal number increase for the consecutive elements, as 


in 32-1, S - 3, 7 - 3S, etc. NMID equals 1 for the node numbering 
increase between successive nodes. 


Free Surface Nodal Cards 
mune ts LO, 4S, 3, Oe, de, 197, 07,96, 86 


Here the total mesh is allowed to move, where the user may fix the 
bottom nodes at any level in the mesh. 


Mid-point Free Surface Nodal Cards 

meee, 1, O8, BS 

Seepage Face Cards 

oe in a 

The seepage is prescribed to start at node 846, NSFC = 84, and 
there are = elements for this seepage face. NSINC is equal to 


as there is a two nodal increase between element corner nodes, 
is 88 —- 86 and 90 - 88. 


bJ 


Hy 
fi 


Statement and End Card 


END OF SILTY SAND DAM. 
eet) 1 


he 





SILTY SAND DAM . 
2,0.01 
m,.,0,6,50,2,1,0 
1,0,0,1 
11,0,110,0 
18, 200,0,1 
28, 200,99,0 
»400,0,1 
, 400, 88,0 
,600,0,1 
§2,600,77,0 
69 ,800,0,1 
79 ,800,66,0 
86,1000,0,1 
96,1000,55,0 
fez, 18,19, 20,13,3,2,4,0.002,0.002,0 
1@,29,35,36,37,30,20,19,4,0.002,0.002,0 
35, 46,52,53,54,4/7,37,36,4,0.002,0.002,0 
92,63,69,70,71,64,54,53,4,0.002,0.002,0 
59 ,80,86,87,88,81,71,70,4,0.002,0.002,0 
lO, 5,2,1 
§6,55,5,2,1 
11,1,28,18,45,35,62,52,79,69,96, 86 
17,34,51,68,85 
36,2,2 
SND OF SILTY SAND DAM. 
3,0.01 
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SILTY S&N0 OAM 


IMEN = 2 FERR «= .91000 


WNPC = 12 NELEMC = 6 WPBOC =: ° NPFS = 6 ITGIY = 59 NPBOC = 2 NSPFAC « 1 ITIME 


GENERATEO NOOAL POINT OATA ANO (X}-(¥) COOROINATES 


NOOE 1 Xz 0.0000 Ys 0.9000 
NOOE 3 Xe ©.90000 Y= 22.9000 
NOOE 6 ) ©. 0000 Y? 44.0000 
NOOE 7 Me ©.0000 Ys 66.9000 
NOOE 9 Xe 0.0000 Ye 68.0000 
NOOE 11 Xe ©.0000 Ys 110.0000 
NOOE 18 Xt 200.0000 Ye 0.90000 
NOOE 20 Ms 200.0000 Ys 19.8000 
NOOE 22 Xt 200.0000 Y= 39.6000 
NOOE 24 Ms 200.0000 Ye 59.4000 
NOOE 26 Xe 200.0000 Ye 79.2000 
NOOE 26 Xt 200.0000 Ys 99.0000 
NOOE 35 Xe 400.0000 Ye 0.9000 
/ NOOE 37 Xs 400.0000 Ys 17.6000 
NOOE 36 Xz 400.0000 Ys 36.2000 
NOOE 41 x8 400.0000 Y= 62.8000 
NOOE 43 Xs 400.0000 Ye 70.4000 
" NOOE 46 Xe 400.0000 Ys 66.9000 
! NOOE 62 Xe 600.0000 Y: 0.0000 
NOOE 54 Mx 600.0000 Ye 15.4000 
NOOE 658 X= 600.0000 Ys 30.6000 
NOOE 56 Me 600.0000 Ys 46.2000 
NOOE 60 Ms 600.0000 Ys 61.6000 
NOOE 62 Ks 600.0000 Y= 77.9000 
NOOE 69 x2 600.0000 Y: 0.0000 
NOOE 71 X= 600.0000 Ye 13.2000 
NOOE 73 Ms 600.0000 Ye 26.4000 
NOOE 75 Xs 800.0000 Ys 393.6000 
NOOE 77 Xe 800.0000 Ye 52.6000 
NOOE 79 Xt 600.0000 Ye 66.0000 
NOOE 386 Xt 1000.0000 Ye 0.9000 
NOOE 66 Xz 1000.0000 Ye 11.0000 
NOOE 90 Xz 1000.0000 Ys 22.0000 
NOOE 62 Xz 1000.0000 Y= 33.0000 
NOOE 64 Xz 1000.0000 Ys 44.0000 
NOOE 66 Xz 1000.0000 Ys 55.9000 
i 
| 
| 
| 
! 
| 
NUMBER OF NOOAL POINTS... 986 
NUMBER OF ELEMENTS....... 28 
ELE. NO. NOOAL POINTS 
) 1 1 12 16 16 20 13 3 2 . 20000E-02 . 20000E-02 °. 
2 3 13 20 21 22 14 6 4 . 20000E-062 . 20000E-02 °. 
3 5 14 22 23 24 16 7 5 . 200008 -02 . 200008 -02 °. 
a 7 18 24 258 28 16 9 6 . 200008 -02 . 20000E - 02 °. 
5 9 16 26 27 268 17 11 10 . 200008 -02 . 200008 -02 °. 
5 16 26 36 36 37 30 20 19 . 20000E-02 . 20O000E-062 °. 
7 20 30 37 36 398 31 22 21 . 200008 -02 . 20000E-02 °. 
6 22 31 39 40 ai 32 24 23 . 200008 -02 . 200008-02 °. 
9 24 32 41 42 43 33 26 258 . 200008 -02 . 20000E-02 ©. 
10 26 33 43 a4 as 34 26 27 . 200008 -02 . 2ZOO00E-02 °. 
11 35 46 $2 63 64 47 37 39 . 200008 -02 . 20000E- 02 °. 
12 37 47 64 56 598 46 39 36 .200008-02 . 2O000E -02 °. 
13 39 46 598 $7 38 49 a1 40 . 20000E-02 . 20000E -02 °. 
14 41 49 58 $9 60 6° 43 42 . 20000E-02 . 20000E-02 3. 
15 43 50 90 51 62 61 46 44 . 20000E-02 . 2O000E- 02 °. 
168 52 53 58 70° 71 94 64 63 . 20000E-02 . 20000E-02 °. 
17 54 94 71 72 73 66 Ce 56 .20000E-02 . 20000E-062 3. 
18 68 9s 73 74 78 96 S65 57 . 200008 -02 . 2O000E - 02 °. 
19 56 66 78 79 77 67 90 59 . 200008 -02 . 20000E -02 °. 
20 $0 67 77 76 78 98 $2 61 - 20000€ -062 . 20000E-62 °. 
21 66 60 6 67 66 61 71 70 . 200008 -02 . 20000E-02 °. 
22 71 61 99 66 60 62 73 72 . 20000E-02 . 20000E-02 °. 
23 73 62 $0 91 92 63 75 74 . 20000€E-62 - 20000E-902 °. 
24 78 53 92 93 64 64 a7 76 . 20000E-02 . 2Z0000E-02 °. 
25 77 64 94 85 98 66 79 76 . 20000E -02 . 20000E-02 °. 


OIRICHLET BOUNDARY CONOITION OATA 


es RESPECTIVE OIRICHLET BOUNOGRY CONDITIONS 
_— 2 110.0000 110.0000 
3- «4 110.0000 110.0000 
S- 6 110.0000 110.0000 
7 6 110.0000 110.0000 
§- 10 110.0000 110.0000 
1- 10 110.0000 110.0000 
§- 687 55 0000 56.0000 
8- $9 $5 0000 56.0000 
d- 81 $5 .0000 55.0000 
%- $3 5B .0000 65.0000 
4- 66 $5 .0000 56.0000 
§- 66 56.0000 66.0000 


DAL POINTS ON FREE SURFACE 
1 17 26 34 ac $1 92 6s 79 66 66 © 


 SEEPAGE FACE NOOAL POINT OATA 
INSEQOUTIYE NOOE NUMBERS ON THE SEEPAGE FACE 


i- 66 | A-18 
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28 
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BOO. 
4oo., 
200. 

oO. 


wry 


oo00 
oo0°0 
oo0°0 
oo00 


<<< <- 
oo as as ca 4 


CONVERGENCE CHECK ERROR §& 


xX COOROINATE 


e900 0 00000 00 


9000 
900 


200 


.0000 
.0000 
.0000 
.0000 
.9000 
.0000 
.00o00 
.9000 
.0000 
. 0000 
.9000 
.0000 
.0000 
.9O000 
100. 
100. 
100. 
200. 
200. 
200. 


200. 
200. 


200 
200 
200 


200. 
200. 
300. 
300. 
300. 
300. 
300. 
300. 
4aoo. 
.0000 
4oo. 
4oo. 
4aoo. 
aoo. 
4oo. 
4oo. 
4aoo. 
4oo. 
aoc. 
soo. 
80°. 
80°. 
80°. 


aoo0 


800 


800. 
80°. 


800 


800 


800. 


800 


80°. 
600. 
800. 
800. 
soo. 
800. 
700. 
.0000 
700. 
700. 
.0000 
7TO0O0. 
BOO. 
8Oo°. 
800. 
8O°. 
80°. 
- 0000 
8O°. 
8O°. 
800. 
8O0°0. 
goo. 


700 


700 


8O°0 


900 
goo 


$00 
9090 
1000 


1000. 


oo00 
e000 
0000 
ooo00 
oo000 
oo00 
oo0°0 
oo0°0 
0000 
oo00 
oo0°0 
oo00 
ooo0o0 
oo0o0 
oo00 
ooo0o 
oo000 
oo0o0 
oo00 
0000 
oo00 


ooo0°o0 
0000 
oo00 
0000 
oo00 
oo0°0 
oo0c0 
oo00 
oo00 
ooo0o0 
0000 
oo0°0 
0000 
ooo00 
oo0°o 
0000 
eRe ze ze) 


oo00 
co00 
0000 
0000 
oo00 
0000 
oo0°0 
0000 
oo000 
0000 


0000 
oo0o0 


oo0o 
e000 
0000 
ooo0°o0 
oo00 
ooo00 


0000 
oo00 
oo00 
ooo00 
o000 
ooo0°o0 


.0000 
goo. 
900. 


1000. 
1000. 
19000, 


1000 


1000, 


1000, 
1000. 
1000. 


1000 


8&8 


oo0°0 
0000 
ooo°o 
ooco 
0000 
oo00 
ooo0°0 
oo0o0 
oo00 
0000 
0000 
ooo0o0 
0000 
oo0°0 


.0c0o0 


¥ COOROINATE 


TOTAL DISCHARGE 


30 


TOTAL DISCHARGE 


rs wee FT 
81.7788 
92.3734 
101.7151 
110.0000 


-1710E-02 


V-POTENTIAL 


END OF SILTY SAND DAM. 


s 


0.0000 110.0000 
11.0000 110.0000 
22.0000 110.0000 
33.0000 110,0000 
44.0000 110.0000 
SS .0000 110.0000 
86.0000 110.0000 
77,0000 110.0000 
88.0000 110.0000 
99.0000 110,0000 

110.0000 110.0000 

©,0000 105.7709 
21.1715 105.7743 
42.3430 105.7846 
63.8145 105.8024 
84.8860 105.8283 

105.8575 105.8638 

0.90000 101.4005 
10.1718 101.4014 
20.3430 101.4044 
30.5145 101.4092 
40.6880 101.4162 
50.8575 101.4248 
61.0291 101.4356 
71.2006 101.4480 
Si3 721 101.4622 
91.5436 101.4780 

101.7151 101.4954 

©.0000 96.8349 
19.4088 96.8392 
38.8177 96.8518 
88.2285 96.8730 
77.8354 96.9026 
97.0442 96.9408 

©.0000 92.0447 

$2273 92.0459 
18.4747 92.0492 
27.7120 92.0549 
38.9494 92.0629 
46.1887 92.0731 
85.4240 92.0855 
84.8614 92.1002 
73.8987 025-3 172 
83.1381 92.1382 
92.3734 92.1876 

0.0000 86.9874 
17.4180 88.9922 
34.8300 87.0086 
52.2451 87.0308 
69.8801 87.0842 
87.0751 87.1073 

0.0000 81.8118 

8.1777 81.8128 
18.3884 81.8188 
24.8330 81.8234 
327 7107 81.8325 
40.8884 81.8441 
49.0891 81.8884 
57.2438 81.8751 
68.4214 81.8944 
73.8991 81.7180 
81.7768 81.7401 

0.0000 75.8434 
18.1314 75.8489 
30.2628 75.8655 
45.3942 75.8932 
80.5256 75.9322 
75.8570 75.9827 

©. 0000 89.6922 

6.9537 89.5935 
13.9074 89.5980 
20.8611 89.6049 
27.8149 69.6151 
34.7886 69.6275 
A1A9223. 99.8431 
48.6760 69.8608 
$8.6297 69.8810 
62.5834 89.7030 
99.5371 69.7274 

©.0000 62.9858 
12.4837 62.8925 
24.9074 62.7129 
37.3811 82.7475 
49.8149 62.7972 
62.2686 62.8634 

0.0000 55.0000 

5.5000 $5 .0000 
11.0000 55.0000 
16.5000 55.0000 
22.0000 55.0000 
27.8000 55.0000 
33.0000 SS .0000 
38.5000 55.0000 
44.0000 55.0000 
49.5000 SS .0000 
S5.0000 SS .0000 

AVE. VEL .161808-03 

: .177798-02 

AVE. VEL .16202E-03 

2 .355988-02 

FERR « .91000 


os 


DOS 


Xe 


a Ve ee wel ce 
ee, ee, ho 


VELOCITY 


~83188-04 
-831S8-084 
83088 -04 
~- 82968-0604 
- 827938 -04 
~§25S82-04 
-82278-04 
-8191E-04 
-81498-04 
.80399F-04 
- 80458-0004 
- 88008 -04 
-8S593SE-04 
-8S$82E-04 
~-8S908-04 
~-8&S30B-04 
-84918-04 
-889SE-04 
.889S282-04 
.889S8-04 
-8895E-08 
.889S8-04 
- 88968-0848 
.8896E-04 
.88988&-04 
.-88998-04 
. 89088-0384 
.89304E-08 
.F3S6E-04 
~-S3S4B-04 
.F3ISOE-04 
-33438-O04 
.833328E-04 
-93318E-O4 
-88O1E°04 
-3938018-04 
~9B8OOE-04 
.9798E-084 
~3I79SE-04 
-8792E-04 
.8378B8E-04 
-37BIE-O84 
-3ITTB8E-O4 
S77Z2E-04 
-37TISE-O4 
~ 1043E-03 
.- 10438 -°03 
~1043E-03 
. 10422 -03 
. 104068 -03 
. 10398-03 
.11088-03 
- 11098-0903 


108E-03 
1088-03 
1088-903 
1088-03 
1048-03 
1038°03 
102E°03 
1018-03 


- 11008°03 
. 12028-0903 
-12028°03 
» 12018-6063 
. 12008-0603 
-11998&-03 
. 113926" 03 
, 1A01E - 63 
13018 °03 


13008 -03 


. 12998 -03 
. 12988-0903 
-12968-03 
.12948-03 
.-1292E°03 
.-12888-03 
- 12868-0903 
-12828-03 
. 14898-0603 
. 14608-0603 
. 14608-0603 
.14628°03 
.1464E-03 
. 14968-0603 
.16188-03 
- 16168-0603 
- 16178-0603 
.16208-03 
.1623E-903 
~1629E-903 
- 1634E-03 
.1642E-03 
. 1649E-03 
.1661E°-03 
-1698E-03 


. 110008402 


.- 1 1QOO00E+O02 


Y-VELOCITY 


112SE-06 


~loe@e = 15 
~-82S18°07 


.17S9E-0O6 


. 29208-0606 


.44998-06 
~-10S88@-15 
.7677E-0O9 
.21339E-67 
.-SS74E-0O6 
~'131&-OS 
-17398-0$ 
~- 24038-05 
~32238-05 
-11868-07 
~-3873E-O6 
~.81678-O6 
_V1S78-05 
- 1630E-05 
. 19098-0O5 
~- 24318-0585 
- 26228-0585 
. 32068-0605 
Ja 295005 
-393648-05 
.6799E-09 
. 87358-0006 
.174SE-05 
.26118-05 
~.3468E-05 
-43318-05 
. 403 7E-08 
.-49218-08 
-93788E-O6 
. 14788-0608 
.19S7E-0O5 
.24S5468-05 
20356 -0Os 
.- 34238-0605 
.- 39086-0905 
-4389E-05 
.- 48748-0585 
.21468-08 
~11138-05 
.22238-05 
.33288-05 
.- 44208-0608 
~-S8O8E-08 
.6S5S948-08 
~-8378E-O6 


_12708-08 
- 19108-0608 
.2S3938-08 
-31718°08 
. 38O2E-08 
.44088-08 
.-SOS3E°-05 
.-§$8898-O0S 
~-6293SE-05 
-32178-09 
.-14618-05 
. 29188-0685 
~4358E-05 
(STEZE-O5 
- 72338-9058 
-161SE-07 


82778E-O6 


-1649E-O$ 
-24618-05 
~3671E-05 
-802468-05§ 
Sasa te-os 
.S4S68-05 


T7O83E-O5 


-96TOB*OS 
.86708-05 
.11O08E-O6 
~17238-08 
. 34768-0058 
.S29308-0S 


7198E-0S 


94378-0858 
-1787E-O6 
-1OS8E-15 
.§$0218-06 
.-10S88E-15 
-1O34E-05 


~1627E-0§ 


. 23198-9058 


~3326E-05 


DO 3 


DO 3 


A-19 


TOTAL VELD. 


~.83188-04 
~-8631S8-04 
~83088-04 
-829368-04 
~-8279E-04 
-825S8-04 
-82278-04 
.- 819 1E-04 
.81498E-04 
~-8O099E-04 
~804SE-04 
.8900E-04 
.8S9SE-04 
~-85838-04 
-8S628-04 
»-8S33E-04 
~-8a497E-084 
.889S8F-04 
.889SB-04 
-889SE-04 
.889368-04 
-8896E-084 
-8898E-04 
. 88998 -04 
-893028-04 
-89058-04 
. 89108-0484 
~-89138-04 
-S$3S6E-04 
.93S55S8B-04 
-933ISZE-O04 
-9346E-04 
-9¥3398-04 
-33288E-04 
~SBO1E-04 
-F8O01E-04 
.8800E-04 
~-3799E-04 
-S797E-04 
.3798E-04 
-37928-04 
-3S7TESE-04 
~978B8E-O4 
-83782E-04 
37778-0848 
- 1043E-03 
- 10438-0903 
.104368-03 
. 10428-0603 
. 10418-0603 
- 10406-0603 
-1108E-03 
. 11088-0603 


- 11O06E-903 
-1108E-03 
.- 11088-0903 
~- 11058-0383 
. 110SE-03 
- 11048-0903 
7 tiga =O3 
-11038-03 
. 11026,*03 
.- 12028-0603 
. 12028-0603 
-1201E-03 
12018-0603 
- 1200E-03 
~- 11998 -O3 
- 13018-0603 
~- 13018 -03 
- 1TJOOE-O3 
~- 12998-0603 
- 12988-0323 
- 12978-0603 
~129SB-03 
~-1293E-03 
.- 1290E-03 
-12878-03 
-128SE-03 
- 14898-0603 
. 14608-03 
. 14618-0603 
_ 14638-03 
.1496E-03 


14698-0603 


-1618E°-03 
- 18188-0603 
- 1617E+903 
- 1620E-03 
-1623E-03 


1629E-03 


~-16348E-03 
-1642E-03 


16S0B-903 


~-1661E-03 
.-1688E-0O3 


-17779E-02 


»-17822E-02 


PRES. 


COEF. 


¢? 





AFPENDIX E-4 


SAMFLE INFUT 


AND OUTFUT 








RUN ADS2 Sra] 2 rer | DEE ies O PEE T ik. 


2,0.0001 

73, 30,0, 30, 20,6,5,0 
1,-340,0,1 
5,-340,34,1 
11,-340,78.0,0 
18,-200,0,1 
22,-200,34,1 
28 ,-200,78.0,0 
35,-80,0,1 
39,-80,34,1 
45,-80,77.0,0 
52,0,0,1 
56,0,34,1 
52,0,74.0,0 
59,50,0,1 
73,50,34,1 

79 ,50,71.0,0 
36,58.75,0,1 
90,58.75,34,1 
96, 59.0,58.0,0 
103,62.0,0,1 
107,62.0,34,1 
113,61.0,58.0,0 


opr oO f- 

< —_e — CE — 
Or or - 
iv vy Oo 
© © 


AWOoOm”mWORDWOO f- 
-— e¢ 


« « « w s se ~~ @ 


Co f~ 





' © — 
oO 

ee, 

rc 
















24,114.0,0,1 
128,114.0,34,1 _ 
134,113.0,60.0,0 
141,118.0,0,1 
145,118.0,34,1 
551,117.0,60.0,0 
$458,122.0,0,1 
}62,122.0,34,1 
168,123.0,60.0,0 
175,126.0,0,1 
179 ,126.0,34,1 
1485,127.0,60.0,0 
192,130.0,0,1 
96,130.0,34,1 
102,131.0,58.0,0 
109,134.0,0,1 
113,134.0,34,1 
119,133.0,58.0,0 
126,138.0,0,1 
30,138.0,34,1 
136,137.0,60.0,0 
'43,142.0,0,1 
47,142.0,34,1 
53,141.0,60.0,0 
60,146.0,0,1 
64,146.0,34,1 
»70,147.0,60.0,0 
77,150.0,0,1 
81,150.0,34,1 
87,151.0,60.0,0 
94,154.0,0,1 
98,154.0,34,1 | 
04,155.0,58.0,0 
m1,156.0,0,1 
15,156.0,34,1 
21,156.0,58.0,0 


-12,18,19,20,13,3,2,4,0.2,0.2,0.0 



















8, me, 35,36, 37,30,20,19,4,0.2,0.2,0. 
5 46,52,53,54,47,37, 36,4, F202, 0. 
2,63,69,70,71,64,54,53,4,0.2,0.2,0. 
9,80,86,87,88,81,71,70,4,0.2,0.2,0. 
6,97,103,104,105,98,88, 87, 4, 0. 20. 


03, m4. 120,121,122,115, 105, 104,4,0. 
20,131,137,138,139,132,122,121,4,0. 
37,148,154,155,156,149,139,138,4,0. 
54,165,171,172,173,166,156,155,4,0. 
71,182,188,189,190,183,173,172,4,0. 
88,199,205, 206,207,200,190,189,4,0. 
05, 216,222,223, 224,217,207,206,4,0. 
22, 233,239,240, 241,234,224,223,4,0. 
39 , 250,256, 257, 258, 251,241,240,4,0. 
56, 267,273,274,275, 268, 258,257,4,0. 
73, 284,290, 291,292,285,275,274,4,0. 
90,301,307,308,309,302,292,291,4,0. 
07,318 ,324,325,326, 319, 309, 308,4,0. 
24,335,341, 342, 343, 336,326, 325,4,0. 
1,352,358,359, 360, 353, 343,342,4,0. 
8, 369, 375,376, 377,370, 360, 359,4,0. 
75,386,392, 393, 394, 387,377,376,4,0. 
92,403, 409, 410,411, 404,394, 393,4,0. 


C000 


0.0 


e 


NN NONNHNONNONNNN NNN NN WN WN WD 
=o oo Qooooo0o0q°n”co0oooqo°o oO 
NNNNONNNNNONN ND NN WN WN LO 


weeeo0o0ooeooeeeoeoeoeooeo°e 


ooooo0oo°oo°o°o°o°oo oO 000 © 


l 
No 






























409, 420,426,427,428,421,411,410,4 
426,437 ,443,444,445, 438 ,428,427,4 
443,454, 460,461, 462,455,445 ,444,4 
(460,471,477,478,479,472, 462, 461,4 
(477,488, 494,495,496,489,479,478,4 
494,505,511,512,513,506, 496, 495,4 
me .0,5,2,1 

96,58.0,1,17,6 

198,58.0,1,17,6 

go0,58.0,1,17,6 

402,58.0,1,17,6 

504,58.0,1,17,6 

eee, 22,45,39,62,56,79,73,96,90,113,107,130,124 

147, fa, 164, 158, Von, 17 Sees, 192, 2159209, 232,226,249, 243, 266, 260 
283,277, 300, 294,317, 311,334,328, 351, 345, 368, 362, 385, 329, 402, 396 
419,413, 436,430, 453,447,470, 464, 487, 481 ,504, 498 
17,34,51,68,85,102,119,136,153,170,187, 204,221 ,238,255,2/2 

289 , 306,323,340, 357,374,391,408,425, 442,459,476, 493 
mo,1,17,198,1,17,300,1,17,402,1,17,504,1,17 

mee o..K = .2..D = 20 FT...END OF PROBLEM. 

3,0.0001 


oo 00 0 ©} 








“se © © &® & - . - ee © 7“ es © © © & & @ . « © &® &® *& -_ e« © &©& &® @& -_ cf e& * . . -_ 2. &® &®& & & @ - = = = - . & * = - 


NOIMEN 


NNPC 


63 


RUN ADS2 


2 


NELEMC 


5P 


390 


e 25° FEET DEPTH 
FERR = »-90010 
NPSOC « ° NPFS 


GENERATED NODAL PDINT DATA AND 


? 
? 
’ 
* 
° 
¢ 
* 
* 
* 
. 
e 
* 
* 
o 
o 
o 
° 
o 
o 
é 
a 
? 
’ 
’ 
’ 
* 
o 
s 
s 
? 
s 
’ 
’ 
’ 
s 
’ 
s 
’ 
? 
. 
? 
’ 
# 
o 
° 
’ 
’ 
’ 
° 
o 
* 
* 
’ 
. 
* 
s 
? 
? 
o 
4 
e 
o 
? 


NODE 

NODE 

NODE 

NODE 

NODE 

NODE 1 
NODE 1 
NODE 2 
NODE 2 
NODE 2 
NODE 2 
NODE 2 
NOD& 3 
NOD& 3 
NODE 3 
NODE 4 
NOoODde 4 
node 4 
NODE § 
NODOE 5 
NODE § 
NODE §56 
NODE 
NODE 6 
NODE 6 
NODE 7 
node 7 
NOOe& 7 
NOOe& 7 
NOD& 7 
NODE 6 
NODE 
NODE 
NODE 9 
NODE 
NODE 5 
NODE 
NODE 
NODE 
NODE 
NOODe 11 
NODE 
NODE 
NODE 
NODE 12 
NODE 
nODe 
nODe 
nODe 
NODE 
NODE 
node 
Node 14 
NODE 
NODE 
NODE 15 


NODE 16 
NODe 
node 
NODE 
NOOe 
NODE 
NODE 17 
NODE 17 
NOOe@ 17 
wODe 16 
NODE 18 
NODE 
NODE 
NODE 18 
NODE 16 
NODE 16 
NODE 20 
NODE 20 
NODE 20 
NODE 21 
NODE 
NODE 
NODE 
NODE 
NODE 
NODE 
NODE 
NODE 
NODE 
NODE 
NODE 
NODE 
NODE 
NODE 
NODE 
NODE 
NODE 
NODE 
NODE 
NODE 
NODE 
NODE 
NODE 
NODE 
NODE 
NODE 
NODE 
NODE 
NODE 
NODE 
NODE 
NODE 
NODE 
NODE 
NODE 
NODE 


16 


17 


22 
22 
22 
22 


23 
23 
24 
24 
24 
24 
24 
25 
25 


25 
26 
26 
27 
27 
27 
27 
26 
26 


29 
29 
26 
29 


31 
31 


160 


2 


164 


1 


173 


6 
7 
6 
1 
6 


160 
182 


4 
& 
6 
Ss 
7 
9 
1 


213 
215 


2 
8 
g 
G 


230 


2 
& 
1 
3 
6 
7 
G 
6 
8 


260 


2 
a 
5 
3 
5s 
7 
d 
1 
3 


260 


2 
4 
5 
6 


300 
307 
306 


1 
3 


315 
317 
324 
326 
325 
330 


xs 
xe 
Xe 
Xe 
Xe 
Xe 
Xe 
Xz 
Xe 
Xe 
Xe 
Xe 
Xe 
Xz 
Xe 
Xe 
Xe 
Xe 
Xe 
Xe 
Xe 
Xe 
Xe 
Xe 
Xs 
Xe 
Xs 
Xe 
Xe 
Xe 
Xe 
Xe 
Xz 
Xz 
Xz 
Xe 
Xz 
Xs 
Xe 
Xs 
Xz 
Xz 
Xe 
Xe 
Xe 
Xe 
Xe 
Xe 
Xe 
Xe 
Xe 
Xe 
Xe 
Xe 
Xe 
Xe 


xe 
Xe 
Xs 
Xs 
xe 
xs 
Xe 
Xs 
Xs 
Xs 
Xe 
Xe 
Xe 
Xe 
Xe 
Xe 
Xe 
Xe 
Xe 
Xe 
Xe 
xs 
Xe 
Xs 
Xs 
xe 
xs 
xs 
Xe 
Xe 
Xe 
Xe 
xe 
Xe 
Xe 
Xe 
xe 
Xe 
Xe 
Xz 
Xs 
Xe 
Xe 
Xe 
Xe 
Xe 
Xe 
Xe 
Xs 
Xs 
xs 
Xe 
Xe 
Xz 
Xe 
Xe 
Xe 
a 
Xe 
Xe 
Xs 
Xz 


- 340 
- 340 
- 340 
-340 
-340 
- 340 
° 200 
- 200 
-200 
- 200 
-200 
-200 
-60 
"690 
“690 


1 1 8 
AAAAA WA ean 
eooogooeo0o e000 00000 


108. 
114. 
114. 
114. 
1133 


Cx) = Y¥) 


.0000 
.-0000 
-0000 
-9000 
.9000 
-9000 
-9000 
. 0000 
. 0000 
0000 
0000 
0000 
- 0000 
. 9000 
. 9000 
-90000 
0000 
-0000 
.0000 
. 9000 
- 0000 
0000 
0000 
. 9000 
. 0000 
.9000 
-9000 
-0000 
.9000 
0000 
- 7800 
. 7800 
. 7800 
.§333 
9167 
.9000 
.9000 
.9000 
-9000 
. $667 
3333 
.9000 
.9000 
.9000 
.9000 
6667 
3333 
.9000 
9000 
0000 
0000 
e667 
3333 
-9000 
.0000 
-0000 


.0000 
~ 3333 
6667 
.0000 
.90000 
- 9000 
.9000 
3333 
. 6667 
-0000 
.0000 
-0000 
.0000 
.3333 
. 6667 
.0000 
-0000 
. 9000 
.9000 
. 6667 
3333 
.90000 
.9000 
.9000 
.9000 
- 5657 
- oe 
9000 
-0000 
0000 
-0000 
- 6667 
. 3333 
-90000 
9000 
-0000 
-0000 
. 3333 
6657 
.9000 
.0000 
-0000 
-0000 
» 3333 
.5657 
.9000 
-0000 
-0000 
. 9000 
. 3333 
- 6657 
.0000 
-0000 
.0000 
.9000 
- 6667 
3333 
0000 
0000 
0000 
0000 
6667 


Vu 
Ye 
Ye 
ye 
Ye 
Ye 
Yu 
Ye 
Ye 
Ye 
Ye 
Ye 
Ye 
Yu 
Ye 
Ye 
Ye 
Ye 
Yu 
Ye 
Ye 
Ye 
Ye 
Ye 
Vu 
Yur 
Ye 
Ye 
ye 
a 
Yu 
ye 
Yu 
Yu 
Ye 
Ye 
Ye 
Ye 
Ye 
Ye 
Ye 
Ye 
Ye 
Ye 
Ye 
Ye 
Ye 
Ye 
Ye 
Ye 
Ye 
Ye 
Yu 
Yu 
Yu 
Yu 


Ye 
Ye 
Ys 
Ye 
Ye 
Ye 
Ye 
Ye 
Ye 
Ye 
Ye 
Ye 
Ye 
Ye 
Ye 
Yu 
Yu 
Yu 
Ye 
Ye 
Ye 
Ye 
Ye 
Ye 
Ye 
Ys 
Yu 
Ye 
Ye 
Ye 
Ye 
Ye 
Ye 
Yu 
Yr 
Ye 
Ye 
Ye 
Ye 
Yu 
Yu 
Ye 
Ye 
Ye 
Ye 
Ye 
Yu 
Yu 
Ye 
Ye 
Ys 
Ye 
Ye 
Ye 
Ye 
Yu 
Ye 
Ye 
Yu 
Ye 
Ye 
Ye 


20 FEET K 


30 


Oo. 
tT. 
34 
48. 
63. 
76 

©. 
17 
34 
4c. 
63 
76. 

°. 
17 
34 
46. 
62. 
TT. 

°. 
17 
34. 
47. 
eo. 
74. 

°. 
17 
34 
46 
$6. 
71 

°. 
17 
34 
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ITGIY 


COORDINATES 


0000 
0000 


.0000 


6557 
3333 


- 9000 


0000 


.9000 
-0000 


6667 


3333 


0000 
06000 
0000 


.0000 


3333 
6667 
0000 
0000 


.9000 


0000 
3333 
6667 
0000 
0000 


.99000 
0000 
~3333 


6667 


0000 


.09000 
.0000 


0000 
0000 


-9000 


0000 


-9000 


0000 


-0000 


0000 


-9000 


0000 
0000 
09000 
6667 


. 3333 


09000 
0000 


.9000 


0000 
e667 


. 3333 


0000 
0000 
0000 


0000 


-6667 


3333 


-0000 
.09000 
.0000 
.0000 
. 5657 
~ 3333 
9000 
.0000 
-9000 
.0000 
.0000 
- 0000 
. 9000 


0000 


.9000 
.9000 
.99000 
.9000 
.99000 
-9000 
-9000 


0000 


. 6667 
Jaa 
.99000 
.9000 
.9000 


0000 


. 6657 
. 3333 
. 09000 
. 9000 


0000 
0000 
6657 
3333 
0000 
0000 
0000 
09000 
5657 


_seaa 


0000 
0000 
0000 
0000 
0000 


.99000 
.0000 
.90090 


0000 


.9000 
. 9000 


90000 


-0000 


0000 
0000 


.0000 
. 6667 


z 


2. 


20 


DEPTH 


NPSDOC 


NSPFAC 


ITIME 


°o 
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: NODE 345 Xs 116.0000 z 34.0000 
: NODE 347 Xs 117.6667 Ys 42.8667 
i NODE 349 Xs 117.3333 Ys §1.3333 
| NOOE 361 38 117.0000 Ys $0 .0000 
| NODE .368 Xs 122.0000 Ys 0.0000 
| NODE 360 Xs 122.0000 ys 17.0000 
L NOOE 362 Xs 122.0000 ys 34.0000 
, NODE 384 Xs 122.3333 ys 42.6887 
[ NODE 3586 Xs 122.6887 Ys 61.3333 
; NODE 368 Xs 123.0000 ys $0.0000 
} NODE 375 Ks 128.0000 Ys 0.0000 
; NDDe 377 Xs 128.0000 ys 17.0000 
; NODe 379 Xs 128.0000 ys 34.0000 
! NDOE 381 Xs 126.3333 Ys 42.6687 
; NODE 383 Xs 128.6667 Ys 61.3333 
, NODE 388 Xs 127.0000 Ys 60.0000 
, NODE 382 Xs 130.0000 Ys 0.9000 
, NODE 384 Ke 130.0000 ys 17.0000 
, NOOE 388 Xs 130.0000 Ys 34.0000 
i, NODE 368 Xs 130.3333 Ys 42.0000 
: NOOE 400 Xs 130.8667 Ys 50.0000 
i NOOE 402 Xs 131.0000 Ys 58.0000 
: NOOE 406 Xs 134.0000 Ys ©.0000 
j NOOE 411 Ks 134.0000 ys 17.0000 
|; NDOE 413 Xs 134.0000 Ys 34.0000 
f NODE 415 Ks 133.8687 Ys 42.0000 
, NODE 417 Xs 133.3333 ys 50.0000 
5 NODE 415 Xs 133.0000 Ys 85.0000 
‘. NODE 428 Xs 138.0000 ys 0.0000 
, NOOE 428 Xe 138.0000 ys 17.0000 
z NOOE 430 Xe 138.0000 Ys 34.0000 
f NOOE 432 Xa 137.6667 ys 42.6667 
i NODE 434 Xs 137.3333 Ys $1.3333 
. NODE 4368 Xs 137.0000 vs 60.0000 
. NODE 443 Xs 142.0000 Ys 2.0000 
j NODE 445 Xs 142.0000 Y: 17.0000 
} NODE 447 Xs 142.0000 ys 34.0000 
x NOOE 446 Xs 141.8887 Ys 42.6667 
#- NOOE 451 Xa 141.3333 Ys §1.3333 
E NOOE 483 Xs 141.0000 Ys 80.0000 
; NDOE 4680 Xs 148.0000 Ys 02.0000 
, NOOE 462 Ks 146.0000 Ys 17.0000 
; NOOE 464 Xs 146.0000 vs 34.0000 
, NDOE 466 Xs 146.3333 Ys 42.6667 
; NODE 488 Ks 148.6587 Ys 61.3333 
1, NODE 470 Kx 147.0000 Ys 60.0000 
.: NODE 477 Xs 150.0000 Ys ©.0000 
., NODE 475 Ke 150.0000 Ys 17.0000 
, NODE 481 Xe 150.0000 Ys 34.0000 
I NODE 483 Xs 150.3333 Ys 42.6667 
: NODE 485 Ks 150.6667 ys 61.3333 
; NODE 487 Kn 151.0000 Ys 60.0000 
; NODE 464 Ks 154.0000 Ys 02.0000 
; NODE 466 Xs 154.0000 Ys 17.0000 
; NODE 465 Xs 154.0000 Ys 34.0000 
; NODE 800 Xs 184.3333 Ys 42.0000 
, NDDE 602 Xs 184.8667 Yr 50.0000 
! NODE $04 Xs 185.0000 ye 58.0000 
, NODE $11 Xs 186.0000 Ys 0.0000 
, NODE $13 Xs 165.0000 Ys 17.0000 
} NODE 6818 Xs 156.0000 Ys 34.0000 
: NOOE 817 Xs 158.0000 Ys 42.0000 
, NODE 8158 Ks 158.0000 Ys $0.0000 
L NODE §21 . 158.0000 Ys $8 .0000 
1 NUMBER DF NODAL POINTS... 821 
NUMBER DF @LEMENTS....... 180 
€Le. NO. NDDAL POINTS 
1 { 12 18 16 20 13 3 2 . 20000E+00 .20000E+00 °. 
2 3 13 20 21 22 14 6 4 . 200008400 .20000E+00 °. 
| 3 5 14 22 23 24 15 7 a . 200008+00 .20000E+00 °. 
4 7 16 24 26 26 18 8 8 . 200008400 .200008+00 °. 
8 8 18 28 27 28 17 11 10 . 200008400 .200008+00 ° 
S 18 28 38 38 37 30 20 18 . 200008 +00 . 200008+00 ° 
7 20 30 37 38 36 31 22 21 . 200008 +00 . 20000800 ° 
8 22 31 38 ao a} 32 24 23 . 200008 +00 .200008*00 °. 
s 24 32 at 42 43 33 28 26 . 20000E+00 . 20000E+00 °. 
10 28 33 43 a4 as 34 28 27 . 200008400 . 20000E+00 °. 
11 35 a6 62 63 S4 a7 37 38 . 200008400 . 20000E+00 °. 
12 37 a7 S4 65 S6 as 38 38 .200008+00 . 20000E+00 ©. 
13 36 a8 S6 87 6s as al 40 .200008+00 . 20000E+00 °. 
14 ai 4s 68 69 sO 50 43 42 . 200008400 . 200008+00 °. 
15 43 50 80 61 62 61 as a4 . 200008400 . 200008+00 °. 
16 82 $3 66 70 71 64 54 83 . 200008400 200008400 °. 
17 64 64 71 72 73 6s 56 S65 . 200008400 . 200008400 ° 
18 66 ss 73 74 75 S68 SS 87 . 200008+00 . 200008+00 °. 
18 Ss Ss 78 78 77 67 SO 58 . 20000E+00 . 200008+00 °. 
20 6° 87 77 78 76 ss $2 61 . 20000E+00 . 200008+00 °. 
21 66 SO Ss 87 88 81 71 70 . 200008+00 . 20000E+00 °. 
22 71 81 6s 8s $0 82 73 72 . 20000E4+00 . 20000E+00 °. 
23 73 $2 90 91 92 63 76 74 . 200008400 . 20000E+00 °. 
24 75 83 92 83 S4 64 77 76 . 20000E +00 . 200008+00 °. 
25 77 84 S4 95 958 68 79 78 . 20000£ 400 . 200008400 °. 
26 86 87 103 104 108 SS 8s 67 . 200008400 , 200008+00 °. 
27 ss 88 108 108 107 6s 90 89 . 20000E +00 . 200008+00 ° 
28 sO 96 107 108 108 100 92 61 . 200008400 . 200008400 °. 
26 92 100 106 110 411 += «101 S4 63 . 200008400 . 20000E+00 °. 
30 64 101 199 +1112 113 102 a5 95 . 200008400 .200008+00 °. 
31 103 114 120 121 122 #115 108 106 . 200008+00 20000800 °. 
32 108 118 122 123 124 #118 4107 108 . 20000E +00 _ 20000E+00 °. 
33 107 118 124 +128 #$128 117 «+1086 108 . 200008400 . 20000E*00 °. 
| 34 109 117 #+126 4%$127: «128~«©118~«6199T~«110 . 200008400 . 20000E+00 °. 
35 111 #118 $4128 #+%129 «#130 «#60119 «©6193~=«112 . 20000E+00 . 20000E+00 oO. 
38 120 131 +137 «+1138 #+%11386 132 #122 «121 . 200008400 . 200008400 °. 
37 122 132 139 +140 «141 «133 «124 ~=«123 . 200008 +00 . 20000800 °. 
38 124 133 141 #142 #143 «134 «#128 © «128 . 200008 +00 . 200008+00 °. 
36 126 134 143 #1446 #148 «#138 128 127 . 20000E+00 . 200008+00 °. 
40 128 #+%138 «#+1465 «#1466 «147 0«13806«6«130~=«128 . 20000E+00 . 20000E+00 °. 
| 41 137 148 184 188 186 149 139 138 .200008+00 . 200008400 °. 
| 42 139 148 188 187 188 180 141 «4140 . 200008400 .200008+00 °. 
43 141 180 188 188 180 181 143 142 . 20000E +00 .200008+00 ©. 
. 44 143 181 180 1861 162 182 148 «41464 . 200008400 . 200008+00 ©. 
| 45 148 182 182 163 164 153 147 148 . 200008400 . 20000E+00 ° 
, 46 154 165 171 #=172 «173 168 186 188 .200008+00 . 200008400 °. 
; a7 158 166 173 174 #178 «#4167 188 187 . 20000E+00 . 200008+00 °. 
a8 188 167 #178 #+$178 «©6177 «»©«188 «©6160 «6158 . 20000E +00 . 200008400 °. 
a9 160 168 +177 «2178 «21768 «189 «69182 = «161 . 20000E +00 . 20000E+00 °. 
| SO 162 169 178 %4180 181 4170 184 «183 . 200008400 . 20000E+00 °. 
S51 171 182 188 189 190 183 173+ «#172 . 20000E+00 . 20000E+00 ©. 
S52 173 183 180 191 192 184 175 174 . 200008400 20000800 °. 
: $3 #175 184 192 163 194 185 177 «#176 . 200008400 . 20000E+00 °. 
. $4 177 #185 194 195 198 186 $179 «+178 . 200008400 . 20000E+00 °. 
| 55 179 18688 196 197 198 187 «+181 180 . 20000E +00 . 200008400 °. 
} S56 188 189 208 206 207 200 180 188 , 200008+00 . 20000E+00 °. 
J $7 190 200 207 208 206 201 192 «191 , 20000E+00 . 20000E+00 °. 
j S58 192 201 209 210 211 202 #194 «1863 . 20000E+00 . 20000E+00 °. BES 
4 59 194 202 211 +212 +213 #203 #196 «195 . 200008+00 . 20000E+00 °. 





$2 207 
83 208 
94 211 
98 213 
98 222 
97 224 
96 228 
99 228 
7° 2306 
71 238 
72 241 
73 243 
74 248 
78 247 
79 258 
77 2868 
78 230 
79 292 
so 294 
81 273 
82 27§ 
83 277 
g4 279 
$s 261 
8? 280 
87 292 
as 294 
88 298 
9° 298 
91 307 
92 308 
93 311 
i 313 
Oe 318 
88 3264 
97 328 
i 328 
98 330 
100° a22 
101 341 
102 343 
103 34$ 
104 347 
108 348 
109 388 
107 390 
106 382 
108 394 
110 3898 
111 378 
112 377 
113 378 
1164 361 
118 383 
119 392 
TAF 394 
118 398 
118 3968 
1290 400 
121 408 
122 aii 
123 413 
124 41§ 
1268 417 
129 428 
127 426 
12§ 4306 
129 432 
13x96 a34 
131 aay 
132 a4as 
133 er 
1364 aag 
135 481 
1398 490 
137 492 
136 aga 
139 498 
1406 406 
141 477 
142 4798 
143 461 
146 483 
146 ass 
1498 an64 
147 4998 
146 436 
148 goo 
180 8O2 


DIRICHLET SOUMDARY 


Oncces 


le 
3- 


° 
° 
° &- 
° 
° 
° 


J+ 

8- 10 

Tier 16 
» 89-102 
2113-102 
- 1968-206 
218-204 
- 300-308 
pati? - 309 
- 402-408 
- 419-408 
 ,-8 04-8106 
@$21-810 


CS 


218 
2198 
220 
221 
234 
236 
239 
237 
238 
281 
2$2 
283 
284 
268 
298 
298 
270 
271 
272 
288 
288 
287 
288 
2688 
302 
3603 
304 
368 
308 
3198 
320 
321 
322 
323 
338 
337 
336 
338 
340 
383 
384 
388 
3$8 
337 
370 
371 
372 
373 
374 
387 
388 
368 
380 
391 
404 
408 
408 
407 
aos 
421 
422 
423 
424 
4268 
436 
438 


440 
441 

442 
ass 
46s 
487 
4$6& 
4s9 

472 
473 
a74 
478 

478 
468 
480 
491 

482 
493 
sos 
807 
sos 
sos 
810 


208 
211 
213 
218 
224 
2298 
2258 
230 
232 
241 
243 
24 
247 
248 
286 
230 
292 
294 
288 
278 
277 
278 
261 
283 
292 
294 
288 
296 
300 
308 
311 
313 
318 
317 
328 
326 
330 
332 
334 
343 
348 
347 
3438 
361 
330 
382 
384 
398 
3386 
377 
378 
38 1 
383 
385 
394 
398 
3968 
400 
402 
ait 
413 
418 
417 
418 
4268 
430 


432 
4364 
4338 
448 
447 
44% 
461 
483 
482 
404 
498 
ass 
470 
479 
461 
483 
466 
487 
4s8 
ass 
$00 
$02 
8O4 


2068 
210 
212 
214 
223 
228 
227 
225 
231 

240 
242 
244 
248 
248 
287 
289 

291 

293 
296 
274 
278 
278 
280 
282 
281 

293 
298 
297 
298 

30§ 
310 
312 
314 
318 
325 
327 
328 

331 

333 
342 
3464 
3438 
346 
380 
388 
391 

383 
33968 
397 

378 
378 
3$0 
3$2 
384 
383 
395 
397 
388 
aol 

a106 
412 
aia 
ais 
ais 
427 
428 


431 
433 
436 
446 
44s 
446 
460 
482 
491 
43 
406 
497 
499 
4768 
480 
482 
ag64 
468 
ass 
487 
488 
801 
803 


RESPECTIVE CIRICHLET SOUNDCARY COWNOCITIONS 


7s .6000 
Ts .9000 
78.0000 
78.0000 
78.0000 
78.6000 
8&8 . 0000 
88.9000 
8$ .coo000 
$8. 0000 
$8 0000 
88.0000 
88.0000 
8&8 .0000 
88 0000 
88.9000 


WOCAL POINTS OW FREE SURFACE 


1 17 
161 187 
361 387 


2 
19 
38 


8 34 
8 204 2 
8 374 3 


217 @2e 228 229 
ois te0 25 Tae 
ays 2268 220 230 
296 290 231 232 
233 239 240 241 
234 241 242 243 
2358 243 244 245 
239 248 249 267 
237 267 2s 269 
280 288 287 288 
251 288 289 290 
282 290 2891 292 
263 202 23 2196 
254 296 298 299 
207 #273 276 «278 
26 2278 ~=«6278 «6277 
7s )60CUd2770=Ci«2780O298 
230 278 280 281 
271 ws: m2 mes 
i644 200 Whi 7392 
288 292 293 294 
209 296 26 29% 
267 299 27 2086 
288 298 299 300 
301 307 308 309 
302 308 310 311 
SOS 311 «+3312 313 
3es 8313 314 318 
306 318 $319 317 
Sos; 32@ 328 320 
318 323 82) sm 
s20 338 338 320 
32% 330 339 sa8 
S22 332 383 334 
235 34m 302 300 
339 343 346 345 
337 348 348 347 
338 347 348 348 
339 349 380 381 
382 388 388 380 
363 3200 381 392 
384 392 393 3964 
39S 394 395 398 
389 398 397 398 
399 #378 $378 377 
3a0 «602377:«63980~= 8 
37% 379 380 361 
372 381 382 383 
373 383 3864 385 
389 392 393 394 
387 396 398 398 
366 309 397 Jes 
388 «630606 «6309 «©6100 
390 400 601 602 
403 408 410 411 
eee 4691 402 613 
GG 413 616 615 
09 618 4189 4617 
an 617 418 ane 
Mee 429 «4627 «3«4%8 
soi 626 428 430 
422 630 631 #32 
423 «82 «35 «ge 
o24 e368 638 46239 
437 463 466 6448 
436 4a6 448 447 
639 667 4668 449 
440 443 480 461 
41 681 682 483 
486 4690 4691 692 
488 492 4693 496 
489 4696 495 «498 
487 499 697 495 
485 498 4699 470 
71 497 e798 G39 
a72 479 «8 6p) 
673 681 482 483 
476 683 484 485 
47S 688 489 4687 
468 694 498 498 
489 46989 497 498 
490 4698 499 $00 
491 $00 $01 6802 
492 802 $03 $064 
$08 811 812 813 
$08 813 816 $15 
SO7 818 819 817 
See 817 B18 5819 
soe $18 $20 6821 
CONOITION OATA 
78.0000 
78.0000 
78.6000 
78.0000 
78.0000 
78.0000 
88.0000 
58.0000 
88.0000 
88.0000 
8%. 0000 
$8. 0000 
$8. 0000 
88.0000 
8S 9000 
88.0000 
as $1 82 9s 
% 228 232 298 
88 3891 402 408 


SEEPACE FACE WOCAL POINT CATA 
CONSEOUTIVE WNOCE NUMBERS CK THE SEEPAGE FACE 


, 88-313 
-198-215 
300-317 
- 402-418 
- 804-821 





ITER 


78 
248 
418 


8s 
288 
4258 


$8 
298 
438 


102 
272 
442 


. Z0000E +00 
. Z20000E +60 
- Z0000E +00 
. Z0000E +00 


s 


ZOOOOE*OO 


. 20000£E +056 
. 200008 +600 
. 200008200 
. 2Z0000£E +00 
. 200008 +00 
. 20000E 60 
. 20000£ 00 
. 20000E 7°00 
. Z20000E 00 
. Z0000£E +60 
. 20000E +00 
. Z20000£E 260 
. Z20000E +00 
. 2ZO000E 00 
- Z20000£ +00 
. ZO0000E +60 
- Z20000£E 00 
. 20000E +00 
- 20000E +60 
. Z20000E +60 
. 290000E 200 
. 20000E +00 
. 20000£E +00 
- 20000E 200 
- 20000E 4606 
- 2Z0000E +00 
- 20000E +00 
- Z0000£ 00 
. 20000£E +60 
. Z0000E +06 
. Z20000E +00 
. Z20000£E 200 
. 20000£ +00 
. Z0000E 200 
. 20000E +60 
. 20000E 00 
. 20000£E 00 
. Z0000E +60 
- 20000E +00 
. Z20000E +00 
- 20000E +00 
- 2O000E +00 
- Z20000£E +00 
. Z20000£E 00 
. Z0000E +00 
. Z0000£E +00 
. Z20000E +60 
. Z20000E 00 
- 2Z0000BE 00 
. 20000£E +00 
. 20000E 00 
. Z20000£ +00 
. 20000E 200 
- ZO000E +00 
- 20000E +00 
- 20000E +60 
. 20000E +00 
. 20000E +00 
. 20000E +00 
. 20000£ +00 
. 2 0000£ +090 


. 2ZOCCOE+ OO 
. 20000E+00 
. ZOO0OE +00 
. 20000E +00 
. 20000E +00 
. ZOCCOCE +00 
. ZOCCCE +00 
. ZOOCCOCE +00 
. 2ZO0O00E +00 
. ZOCOOE +00 
. 20000E +00 
. ZOOO0OE +00 
. ZOCOOE +00 
. ZOCOCOE +00 
. 2OOOOE +00 
. ZOOCOCE +00 
. ZOOOOE +00 
. 2OCOCOE +00 
. 2OCCOE +00 
. 20000E +00 
. 20000E +00 
. 20000E +00 
. 20000E +00 


113 
283 
483 


118 
288 
488 


1390 
300 
470 


2 OOOOE +00 


- 20000E 400 
. 20000E +00 
. Z20000E +00 
. 20000E 00 
200008 400 
. 20000E 200 
. Z0000E 400 
. 20000E +00 
. 20000E +00 
. Z0000E +00 
.20000E 400 
. 20000E +00 
. Z0000E +00 
. 20000E +00 
. 20000£ +00 
. Z0000£ 400 
.20000E400 


+ 


2O000E +00 


. 20000E 200 
. 20000E 400 
.- 20000E +00 
. 20000E +090 
- 20000E +00 
.- 20000E +60 
.- 20000E +60 
.- 2ZOGCOCOCE +00 
- 20000E +090 
- Z0000E +00 
. 20000E +00 
. ZO000E¢909 
- 20000E +090 
- 20000E 00 
. 20000£E +00 
. Z0000£E +00 
- Z0000£E +00 
- Z20000E +60 
- 20000E +60 
. 20000£E +00 
- 20000£E +00 
- 20000E 200 
. Z0000E 200 
. 20000E +00 
- 20000E 00 
. 20000E 00 
. 20000E +00 
. 20000£ +00 
. Z0000E +00 
- 20000E°00 
. Z0000£E +00 
- 20000E +00 
. 20000E +00 
. 20000E 200 
- Z20000£E +00 
- 20000£ 00 
. Z0000E£E +00 
. 29000E 00 
. 2O000E 00 
. Z0000E +00 
- 20000E +00 
. 20000£E +056 
- Z20000£E +00 
. Z0000E 4056 
. Z0000E +656 
- Z20000E +00 
- 20000E 400 


- Z0000E 400 
. Z0000E 200 
. Z0000£ +00 
. Z0000E +00 
. ZO000E +00 
. Z0000E 200 
. Z0000E 400 
. 20000£E 400 
- 20000E 00 
. Z0000E +60 
. Z0000€ +60 
. Z0000£ +00 
. ZOO00E +00 
. Z0000E +00 
. 20000E 200 
. Z0000£ 400 
. 20000E 200 
. Z0000E +00 
. 200002 +00 
. 20000£ +00 
. Z0000E +00 
. Z0000£ +00 
. 20000£ +00 


138 147 
3068 317 
478 487 


183 
323 
493 


000000000000000000000000000000000000000000000000000000000000000000 


eoo0oo0ogo0go0coooooo0o0oco0o0000000090 


1964 
334 
$04 


1790 
340 
© 


B-6 





NOOE«457 xs 180). 6272 Y « $8.1055 
NDDE2470 x = 148.9281 Voue 88.1288 
 ~NDOE483 Xo: 141.0718 Y « 85.1328 
_ NODE?438 x © 137.0728 Ys 88.1120 
+ NDOE418 x 2 133.0000 Ys 58.0000 
NODE2402 xX = 131.0000 Ys 88.0000 
NODE2388 x s 128.8874 Ys 88.1382 
NODE: 388 Xs 122.8288 Y = 88.1999 
NODE:351 x s 118.7978 Y 3 55.3181 
NODE2334 xX © 113.0710 Y 3 55.1845 
NODE?317 x 2 109.0000 Ys 58.0000 
NODE2300 x 8 107.0000 Ys 58.0000 
NODE2283 x 2 102.9313 Y « a8.2148 
NODE2295 x3 98.9338 Y 3 $8.2707° 
NODE2248 Mos 93.0588 Y « 88.2937 
NODE2232 x o3 59.0571 Yo: 55.2580 
NODE2215 x. 85.0000 Ys 58.0000 
NODE:198 Mo 83.0000 Y * 88.0000 
NODE?181 x8 78.9382 Y 3 55.3828 
NODE2164 X38 74.9424 Y « 58.5028 
NODE2147 x oa8 89.0554 Y « aa.56805 
NODE=130 x8 88.0578 Ys 58.4979 
NOOE©113 MF 61.0000 Ya 58.0000 
NODE: $8 Mo 89.0000 Ys 55.0000 
NODE? 78 Kos 50.0000 Ys as.5559 
NODE: 682 x3 0.0000 Y « 52.3818 
NOOE: 48 KX 2 +*80.0000 Y : 58.1508 
NOOE: 28 X © -200.0000 Ys 71.5308 
NODE: 11 X © °©380.0000 Y « 78.0000 


ITER 3 2 





NEW CDOROINATES OF THE FREE SURFACE LINE 
NOOE:487 x 3 180.9255 Ys aa.1412 
NOOE#470 x 3 148.9289 Y 3 88.1778 
NOOE?483 Ks 141.0899 Ys as.1821 
NOOE2436 KM os 137.0712 Ys aa.1489 
NOOE:419 xX os 133.0000 Yo" 88.0000 
NOOE2402 x 8 131.0000 Ys 88.0000 
NODE2368 Ks 125.8023 Ys 88.1849 
NOOE:388 Kos 122.9318 Ys 88.2272 
NODE:381 KX of 119.7315 Ys 88.3038 
NODE2334 x 3 113.0893 Y « 88.1881 
NOOE:317 x 3 108.0000 Y 3 58.0000 
NODE!300 x 3 107.0000 Ys 58.0000 
NODE?283 x 3 102.9334 Y «x 88.2872 
NODE?288 xX 8 88.9355 Ys as.33ai14 
NODE2248 xX 3 93.0922 Ys 88.3838 
NODE?2232 x 3 89.0943 Ys 88.3291 
NODE2218 x 3 88.0000 Yo 88.0000 
NODE2188 x 3 83.0000 Ys 58.0000 
NODE:181 x 3 78.9410 Ys 58.4882 : 
NODE:184 x 3 74.9499 Ys 58.9188 
NODE?2147 x 3 58.0803 Y 3 58.8934 
NODE2130 x 3 68.0538 Y 3 58.9088 
NODE*«113 KX 3 81.0000 Y 3 88.0000 
NOOE: $9 x 3 8&8 .0000 Ya 88.0000 
NOOE: 79 x 3 80.0000 Y 3 88.8918 
NOOE: $2 x 3 0.0000 Y 3 82.8589 
NODE: 48 KX 2 =80.0000 Y 3 99.8218 
NODE: 28 X 2 *200.0000 Y 3 72.1278 
NODE? 11 X 2 *340.0000 Y 3 78.0000 
4 ITER «= 3 
é 
' 
NEW COORDINATES OF THE FREE SURFACE LINE 
NODE:487 xX 3 150.9286 Ys 58.1410 
NODE:470 x 3 148.9286 Ys 58.1772 
NOOE=483 x 3 141.0998 Y « 58.1818 
_ NDOE«4&38 x 3 137.0712 Ys 58.1488 
NODE«418 x 3 133.0000 Y 3 58.0000 
NOOE #402 x 3 131.0000 Y « 58.0000 
NOOE:388 x 3 128.9015 Y 3 88.1834 
NODE:3868 x 3 122.9318 Y 3 $8.2271 
NODE:381 xX 3s 118.7340 Y 3 88.2399 
NODE:334 x 2 113.0880 vi) = 88.2048 
NOOE:317 x 2 109.0000 Y « 88.0000 
. NODE? 300 x 3 107.0000 Ys 88.0000 
NODE?283 x 3 102.8333 Ys 88.2871 
NODE:288 x 8 95.8386 Ys 85.3808 
NODE 249 x 93.0822 Y 3 §8.3822 
NODE1232 Xs 59.0843 Y 3 §5.3250 
NOOE?215 a 55.0000 Y 3 58.0000 
NOOE? 186 Kos 53.0000 Y 3 88.0000 
NODE?#181 x 3 78.9410 Y 3 58.4955 
NOOE:184 x 3 74.9499 Y 3 88.9183 
NOOE?147 x 3 99.0804 Y 3 88.8908 
NOOE2130 x 3 95.0538 Y 3 88.8082 
NODE2113 x 3 51.0000 Y 3 58.0000 
NODE: 898 x 3 59.0000 Y 3 88.0000 
NODE: 78 x 3 80.0000 Y 3 58.5784 
NODE: $2 x 3 0.0000 Yo: 82.5584 
NODE: 45 KX = =80.0000 Yo S§6.5380 
NODE: 25 KX © +-200.0000 Y 72.1238 
NODE: 11 X 2 *+©340.0000 Yo: 78.0000 
CONVERGENCE CHECK ERROR 3 . 82272-0858 
ONDDE X COOROINATE Y COOROINATE V-POTENTIAL X-VELOCITY Y-VELOCITY TOTAL VELO. PRES. COEF. 
1 -340.0000 0.0000 78.0000 -.8214E-02 -. 1203E-04 .8214E-02 
2 -340.0000 8.8000 78.0000 -.8211E-02 - .9783E-14 .8211E-02 
3 -340.0000 17.0000 78.0000 -.8202E-02 .8115E-08 .§202E-02 
4 - 340.0000 25.5000 78.0000 -.8188E-02 }. .8189E-02 
5 - 340.0000 34.0000 76.0000 -.8188E-02 . 2002E-04 .818SE-02 
9 -340.0000 41.3333 78.0000 -.8141E-02 -. B8870E-14 .8141E-02 
7 -340.0000 45.8867 78.0000 -.8112E-02 .3081E-08 .8112E-02 
8 -340.0000 55.0000 78.0000 - ,8077E-02 °. . 8O077TE-02 
9 -340.0000 83.3333 78.0000 =, 8038E-02 .4931EF-04 .8038E-02 
10 -340.0000 70.8887 78.0000 -.7987E-02 .17018-13 .7987E-02 
| -340.0000 78.0000 78.0000 - .7938E-02 . T8B7E-04 .7938E-02 
12 -270.0000 0.90000 78.0770 - .8490E-02 -.2288E-05 .8480E-02 
13 -270.0000 17.0000 78.0801 -.86484E-02 .830SE-04 .848SE-02 
14 -270.0000 34.0000 75.0894 -.8488E-02 .1287E-03 .8488E-02 
18 -270.0000 47.8873 78.1019 -. 8447E-02 .1844E-03 .8448E-02 
18 -270.0000 61.3748 78.1194 -.58418E-02 . 2483E-03 .§84228-02 
17 -270.0000 75.0819 73.1427 -. S353E-02 .3210E-03 . 8388E-02 
18 -200.0000 0.0000 72.0870 -.8784E-02 .1814E-08 .8784E-02 
18 -200.0000 8.5000 72.0879 -.8784E-02 44886-0484 .8784E-02 
20 - 200.0000 17.0000 72.0808 -. 8753E-02 . 8456E-04 .8784E-02 
21 -200.0000 25.5000 72.0888 -.8783E-02 .1343E-03 .8784E-02 
22 -200.0000 34.0000 72.0723 -.8782E-02 .1854E-03 .8754E-02 
23 -200.0000 40.3840 72.0783 - .8783E-02 .2084E-03 .S75SE-02 
24 -200.0000 45.7078 72.0888 -.8784E-02 .2872E-03 .8788E-02 B=] 
25 -200.0000 63.0519 72.0935 -.5788E-02 .2997E-03 .8790E-02 
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93996 
3691 
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© .90000 
17.0000 
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45.7636 
67.8573 
69.3606 

0.0000 

6.6000 
17.0000 
28.5000 
34.0000 
38.4387 
44.8793 
60.3190 
85.7557 
61.1663 
68.6360 

0.0000 
17.0000 
34.0000 
44.2206 
54.4411 
64.6617 

©o.0000 

6.6000 
17.0000 
28.5000 
34.0000 
36.7606 
43.6815 
46.3427 
§3.1236 
67.6046 
62.6854 

0.0000 
17.0000 
34.0000 
43.0601 
62.1203 
61.1504 
©.90000 

§ .6000 
17.0000 
28.8000 
34.0000 
338.2762 
42.5565 
46.8377 
61.1168 
66.3662 
56.6764 

0.0000 
17.0000 
34.0000 
42.2762 
60.8686 
66.8377 
0.0000 
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17.0000 
25.8000 
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80.2021 
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0.90000 
6.69000 
17.0000 
26.6000 
34.0000 
35.1010 
42.2021 
46.3031 
80.4041 
64.6051 
66.6062 
0.90000 
17.0000 
34.0000 
42.2161 
60.4323 
66.6454 
©o0.0000 
6.8000 
17.0000 
26.6000 
34.0000 
38.1181 
42.2302 
46.3463 
60.4604 
64.8786 
65.6906 
0.0000 
17.0000 
34.0000 
42.2182 
$0.41363 
535.8845 
0.0000 
8.5000 
17.9000 
25.5000 
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4136 
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- 4166 
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3164 
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.O761 
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7437 
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2863 
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4487 
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- 6650 


. 4167 
. 3270 
.2185 
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-7241 
- 6082 
- 6746 
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. 3505 
. 1669 
- 6870 
-6710 
- 6120 
.6933 
- 5382 
- 4369 
. 2555 
. 2006 
-O960 
- 33800 
-§831 
. 7367 
. 6607 
-8811 
.4783 
2294 
-O386 
.8088 
-8762 
- 4922 
4736 
4167 
- 31854 
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.- 91668-0602 
.9184E-02 
-91608-02 


-9 1446-02 


-91368-02 
-9126E-02 
.95358&-02 


-6832E-02 
-93829E-02 


95228 -02 
-95148-02 


-S6O07E-02 
-6500E-02 
-9492E-02 


9452E-02 
9476E-02 


-94688-02 
-8763E-02 
-6779E-02 
-9616E-02 
- 86358-0002 


S6S6E-02 


.93556E-02 
- 1004E-01 
- 10068-01 


10118-01 
1018E-01 
- 10238-0601 
10228-01 


- 10188 -01 
~- 10118-0171 
- 10008 -01 
.6934E-02 
-9846E-02 
.§317E-02 
- 95808-0602 
- 1032E-01 
. 1096E-01 
.-11718-01 
. 1200E-01 
- 62462-0602 
-S368E-02 
. 66628-90902 
. SOG6BE-02 
- $85 7E-9O2 
. 10608-01 
- 11438-0601 
- 12688E-01 
. 14608 -01 
- 1676E-01 
-12688-01 
-7T7O8E-02 
- 738418 -02 
- 8862E -02 
-1017E-01 
- 4478-01 
-3732E-01 


TFIOIE-O2 


-T3IT8E-O2 
- 7478 E-02 


T776E-02 
®1638-02 
84666-0602 


- 88238-0602 
-6878E-02 
. 10148-01 


. 17608-01 
_- 31468-0601 
- 71608 -02 
- T314E-02 
- 7648E-02 
.6323E-02 
-6830E-02 
-16628-12 
- 68628-0902 
. 7045 E-02 
- 71228-0902 
.- 73118-°02 
- 7T827E-02 
- 78828-0902 
- 78418-0602 
-76881E-902 
- 68658-0902 
-9328618-03 
-2283E-01 
- 86028 -02 
- 86028 -902 


71&42E-02 


- TOBLE-02 


&462E-02 
31486 -01 
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. 8826E-02 
. 8873E-O2 
-6736E-O2 
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6671E-02 


-6646E-02 


S4356-02 


- 41426 -02 
-37T3ISE-02 
-6764E-02 
. 63678-0602 
.-64368-02 
. €408E-02 
-S967E-02 
- 358528 -02 
.42508-02 
.6162E°-02 
.-81988-02 
- 62188 -02 
-817SE-02 
-SOBBE-02 
-S676E-02 
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-&4&724E-02 
- 33808-0602 
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-S754E-02 
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-2023E-03 
. 27198 -03 
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.41278-03 
-SO36E-05 





S471E-04 


.11968-O3 
. 16648E-03 
-2378E-03 
.- 26018-03 
-31198-03 
. 32608-03 


3520E-03 


- £006£-O3 
.4827E-0O3 
. 89618-0858 
. 809 38E-04 
- 1668E-9O3 
. 2708E-03 
- 36648-0383 
-4763E-903 
- 367OE-04 


22128-0484 


-6914E-05 
- 5 4648-06 
-17208-03 
. 74718-0848 
- 3S6O03E-O3 
. 18668-0603 
.§80088-03 
- 38438 -03 
.6465E-03 
. 17638-04 
- BO04EB -O3 
- 96 83E-03 
-S6TOJE-03 
- 43282-0603 
-54248-03 
- 35 238-04 
. TO76E-O3 
. 13678-0602 
. 22738 -02 
. 30178-02 
- 3SE67TOE-O2 
- 36618E-02 
_. 406 7E-02 
- 33I33E-O2 
- 168868 -02 
- 44568-0902 
- 16478-0484 
- 18808 -02 
-3817E-02 
-51868-02 


* 


* 


8346E-02 
4608E-03 


-6€302£-08 
. §268E-03 


16612-90902 
27638-02 
3688E-02 


-4678E-62 
-§8820E-02 
.8385E-02 


* 


8827E-02 


- 2644E-01 
- 44426E-01 
- 14738-0484 
- 1676E-O2 
- 40328 -02 
. 86438E-02 
-8118E-02 
-4682E-01 
- 1838E-084 
-86638-03 
.17048-02 
- 26646 -02 
- £06 6E-02 
-8144E-02 
-6018E -02 
86628 -02 
- 8683E-02 
- 28888 -01 
-42188-01 
- 13648-0604 
- 1I7FIIE-O2 
-4142E-02 


6O86E-02 
6624E-02 


. 1O508-01 
. 14146-0848 
- 58158 -03 
- 176588-02 
-2821E-02 
-4156E-02 
. SOBBE-02 
. SOB3E-O2 
~-7T331E-02 
- 95823E-02 
- T664E-02 
- 10458-02 
-1264E-04 
~ 17718 °02 
. 41408 -02 
- €6620-02 
-61688&-02 
-11388-02 
. 12648-0848 
. 88158-03 
- 1761E°02 
-2814E-02 
-41088-02 
- 4685 78-02 


S633E-02 


-61318-O2 
. 6679F-9O2 


33890E -02 


- 44968 -03 
- 10698-0004 
-17838-02 
- 40818-0602 
-SIVITE-O2 
-S776E-02 


2478E-03 


- 1064E-04 
.88618F-O3 


17SOE-02 
28618-0902 





-89166E-0O2 
.9185E-062 
.8162E-02 
.91458-02 
.9142E-02 
-9134E-02 
96 36E-02 
- 5338-02 
.S530E-02 
. 8823E-62 
.8617E-02 
. 88108-0602 
- 360SE-02 
- 949 7E-02 


34380E-02 


-9484E-02 
.6479E-02 
. 3 763E-02 
.9778E-02 
.9818E-02 
.9642E-02 
. 96638-02 
-93567E-02 
- 1004E-01 
- 1008E-01 
- 1011E-O1 
-10168-01 
- 10238-0901 
- 10228-01 
_- 1019E-0O1 
_ 1011E-O1 
. 10028-01 
.9940E-02 
. 85866E-02 
.9317E-02 
. 9566E-02 
_ 103 7E-01 
.- 11028-01 
- T172E-01 
-1202E-01 
-6246E-062 
-6368E-02 
-6876E-02 
. 9376E-02 
- 1031E-O1 
-11228-01 
- 1207E-0O1 
.- 13618-01 
- 1617E-01 
. 158878-01 
- 1363E-01 
- T7TOB E-02 
. 808 1E-02 
- 86008-0002 
.11418-0O1 
.18718-01 


-37328-°01. 


TIOIE-O2 


- 7T&24E-02 
- 7TE86E-02 
-§252E-02 
- 8068E-O2 
-8838E-02 
- 1083E-01 
- 12808-0901 
-1344E6E-01 


- 33488-01 
- S443E-01 
- 7T16OE-02 
- TEOLE-O2 
- 588218 -02 
- 10238-01 
-1207E-01 
- 4682E-01 
- 6662E-02 
. TOB TE-O2 
.- TIZIE-O2 
. 788 2E-02 
. 8886E-02 
- 92288-0902 
-8883F-02 
.-11446-01 
-12288-01 
- 28818 -01 
-4787E-01 
. BEO2ZE-02 
-71186E-02 
.6286E-02 
.- 63246-0002 
. 10898-01 
- 33168-0901 
. 869 6E-02 
- 6884E-02 
- B900E-02 
- 7342E-02 
.- 79S 3E-02 
.-6391E-02 
-6923E-02 
.8126E-02 
- 1066E-01 
-8826E-02 
.6646E-02 
. 639 7E-02 
._ 6676E-O2 
- TEJOE-O2 
-SITZE-O2 
-8826E-02 
- 4400E-02 
-61862E-02 
-6262E-02 
-S€486E-02 
-685828E-02 
- 73428-0602 
-7T6E26E-O2 
. 78S 7E-02 
- TT40E-02 
- T477E-02 
- 41366-0602 


e 


7641E8-03 
S9S88E-02 


-6238E-02 
- 7038E-02 
- 7825E-02 
- S9SSE-O2 
- 2472E-02 


S754E-02 
SS39E-02 
6O28E-02 
63IZ27E-O2 


B-8 


we 










161 74.4734 46.3091 68.8643 -. 4663E-602 - SI3S3E-O2 - 7084E-02 
| 162 74.6312 80.4122 $8.7810 -.4516E-02 - .6408E-0602 - TO€&BE-02 
193 74.7830 84.8182 88.6873 -,49189E&-02 -.3046E-02 -8786E-02 
164 74.9469 88.6183 88.6183 - .6217E-02 *. 1124E-03 -8219EF-02 
198 76.0005 ©0.90000 89.4388 -.8878E-02 -9OIJOCE-O8 . 8576E-02 
168 76.9000 17.0000 89.3603 -.6621E-02 -. 17TOE-02 -6768E-02 
167 76.0000 34.0000 89.1203 - .SIT1IE-02 -.3612E-02 - 9484E-02 
168 79.3146 42.1808 66.9293 - .4722E-02 -.8061E-02 -6921E-02 
169 79.6263 60.3613 88.7084 -.4695E-02 -.6771E-02 . TAROE-0O2 
| 170 76.9436 86.5419 86.6863 -.7668E-02 -.4188E-03 - T6E79E-02 
171 78.0000 ©.0000 89.3807 - .63I78E-02 -S494E-05 -S376E-02 
172 78.0000 6.8000 86.3623 - .S8348E-02 -.8764E-03 -6420E-02 
173 78.0000 17.0000 $9 .3062 -.6312E-02 -.1786E-02 . 8894E-02 
174 78.0000 26.8000 69.2101 -.8141E-02 -.2779E-02 -8644E-02 
178 76.0000 34.0000 59.0700 - . &4907E-02 -. 3832E-02 .6226E-02 
176 78.1898 38.0778 88.6824 -, &@690E-062 -.4&381E-02 .6418E-02 
177 78.3137 42.1882 68.8841 -.4481E-02 -.4981E-02 . 8887TE-02 
178 78.4706 46.2327 538.7761 - .48478-02 - .S317E-02 . 6997E-O2 
179 76.6273 60.3103 68.6601 - .4624E-02 -.8238E-02 . T7B63E-O2 
180 78.7842 54.3879 88.8341 -.8267E-02 -.4486E-02 -63G2E-02 
161 78.9410 88.4688 $8 .4984 -.1058S5E-01 -. 28138-0603 . 1088E-01 
182 80.0000 ©.0000 66.3261 -.8168E-02 .S969E-05 .8169E-02 
183 80.0000 17.0000 889.2844 -.S080E-02 -.1736E-02 .836GE-02 
184 80.0000 34.0000 88.0224 - .4637E -02 -.3744E-062 . S880E-02 
188 80.3235 42.0776 68.8421 - .#4231E +02 - .4799E-02 - 639 8E-02 
186 80.6470 §$0.1862 86.6133 -.4919EF-02 -.8923E-02 - TESSE-O2 
187 80.9706 68.2327 88.3186 -.2388E-01 -.8619E-02 -2424E-01 
r 188 82.0000 0.0000 69.2774 -.4977E-02 -6271E-086 49 77E-02 
188 82.0000 8.5000 89.2863 - .4936E-02 -.8888E-03 -BOOSE-02 
190 82.0000 17.0000 86.2046 -.4880E-02 - .1713E-02 -B81TF2E-02 
191 82.0000 28.8000 68.1116 -.4969E-02 -.29878E-02 -S$381E-02 
162 82.0000 34.0000 86.9773 -. &4383E-02 -.38S1E-02 -8712E-02 
193 82.1967 38.0000 68.8983 -.4158E-02 -.#191E-+02 -8801E-02 
194 62.3333 42.0000 68.8027 - .3692E-02 -.4631E-02 - BOS49E-02 
195 $2.8000 48.0000 $8.3984 - . 3800E-02 -.8861E-02 -8818E-02 
198 82.8867 6o.0000 88.8700 -.3148E-02 - .&48S80E-02 -8818E-02 
197 82.8333 54.0000 88.4098 -.86949E-02 -. 1386E-01 - 13SS9E-01 
168 83.0000 §8.0000 88.0000 -.1637E-01 -.2302E-01 -2826E-01 
166 84.0000 0.0000 $9.2267 - .&77T1E-02 . 4780E-08 -&TT1E-02 
200 84.0000 17.0000 59.1870 -.48S56E-02 -. 1687E-02 - 4962E-02 
201 84.0000 34.0000 $8 .9347 -.4140E-02 -.36886E-02 - S45 7E-02 
202 84,0000 42.0000 $8.7715 - .3847E-02 -,@4B27TE-02 -8813E-02 
203 84.0000 80.0000 688.8815 - .2934E-02 - .47998&-02 -8628E-02 
204 84.0000 88.0000 88.0000 - .9OSBE-~13 -.2469E-01 - 2499E-01 
208 83.0000 ©.0000 89.1818 -.48889E-02 - 3847E-O8 -4886E-02 
206 86.0000 8.8000 89.1644 -.4835E-02 -.8266E-03 -4910F-02 
207 86.0000 17.0000 86.1118 - .4468E-02 -.195S8E-02 -4766E-02 
208 $6.0000 28.8000 88.0221 -.4228E-02 -.2683E-02 -4936E-02 
208 66.0000 34.0000 88.8946 -.3921EF-02 -.3488E-02 -8226E-02 
210 68.6333 38.0000 88.8235 -.3888E-02 - .36899E&-02 -8402E-02 
211 66.9687 42.0000 §6.7420 -.3424E-02 -. 44086E-02 -6S7E6E-02 
212 66.6000 43.0000 $6.9480 -.2908E-062 - .S$386E-02 -8132E-02 
213 68.3333 80.0000 66.8308 -.3012E-02 -.4666E-0602 - 8801E-02 
214 66.1697 94.0000 88.3730 -1881E-02 -.1321E-01 . 1IJOE-O1 = 
216 66.0000 88.0000 $8 .0000 -1168E-01 -.2136E-01 - 2449E-01 
216 88.0000 ©o.0000 66.1372 -.4383E-02 - 2346E-08 -&3B6IE-O2 
217 68.0000 17.0000 66.0880 -.4263E-02 -.1926E&-02 -4964E-02 
216 88.0000 34.0000 68.8999 -. 3687E-02 -. 3363E-062 - 4684E-02 
2198 87.9774 42.0847 66.7083 -. 3101E-02 -.4263E-02 .6272E-02 
220 87.3546 60.1083 66.8124 -.1281E-02 -.6824E-02 .8671E-02 
221 87.0322 686.1640 88.2375 .1961E-0O1 -.8849E-02 . 1791E-01 
: 222 90.0000 0.0000 $9 .0943 - .42088-02 . 10898-08 .4208E-02 
i 223 90.0000 8.6000 $8 .0773 -.41818-02 -.79 798-03 .42278-02 
224 90.0000 17.0000 98.0268 -  &OTTE-02 -.18688-02 -4378E8-02 
| 226 90.0000 28.8000 68.9412 -. 36228-0902 -.2420E-02 .4824E8-02 
226 60.0000 34.0000 68.8208 -. 34668-0602 -.3248EF-02 .4774E8-02 
° 
227 89.8441 38.0847 $8.7333 -. 32086 8-02 -,.3692E-02 - 4896E-02 
| 228 $s .6881 42.1083 86.6773 -.2607E-02 -.4088E-02 .$026E-02 
228 86.9322 49.1940 $6.8918 -.21888-02 -. 44828-0602 - 49 78E-02 
230 $9 .3762 80.2187 66.4686 -.13828-02 - .8381E-02 -8$8S6E-02 
231 86.2203 84.2734 86.3616 . 2801E-02 -.4118E-02 - 4618EF-02 
232 88.0643 $8.3280 86.3280 .82108-02 -. 76028-0603 -§$266E-02 
233 92.0000 ©.0000 $89 .0833 - .4008E-02 -.82368-08 - $008E-02 
234 62.0000 17.0000 86.6868 - .3872E-02 -.18638-02 . 41768-0602 
238 $2.90000 34.0000 88.7866 -. 32638-0602 -.31418-02 . 48818-0902 
238 91.9878 42.1184 88.6463 - .2TOOE-02 -.3860E-02 .4738E-02 
237 91.3786 80.2396 88.4888 -.12388-02 -. 48748-0602 .4736E-02 
236 91.0833 88.3881 $8.3978 -27208-02 -.62648-03 .27918-02 
238 94.0000 ©.0000 86.0141 - .38378-02 -. 2Z2307E-08 -38378-02 
240 94.0000 8.8000 88.393678 - .3763E-02 -.7989E-03 - 38GOB-O2 
241 64.0000 17.0000 $8.9490 -.37TO8E-02 -.19308-02 -4011F-02 
242 64.0000 28.8000 $6.6680 -.3461E-02 -.2284E-02 -4136E-02 
243 94.0000 34.0000 $8 .7846 -.3132E-02 - .3037TE-02 -4362E-02 
244 63.6437 38.0837 88.3622 -. 28948-0602 - . 33878-0602 - 44138-0902 
2458 63.6674 42.1274 98.6226 -.2879E-02 -. 39938-0602 -4802E-02 
246 93.8311 46.1911 88.8486 - .2028E-02 - . 37898-0602 . 42718-0602 
247 63.3745 80.2848 88.4734 -.1288E-02 -. 36268-0602 - 40408-0602 
248 63.2186 4.31898 68.4095 -. 1TOOE-0O3 -.2246E-02 .2289E-02 
248 63.0622 $8 .3822 $6 .3822 . 8826E-03 -.23338-03 .- F72Z213E-0O3 
280 66.0000 ©o.0000 $8.93798 - .3648E-02 -.4886E-08 -3648E-02 
281 66.0000 17.0000 $8.9130 -.3818EF-02 -.14689E-02 .- 38218-02 
282 96.0000 34.0000 $8.7246 -.2¥9S8E-02 -.28936E-02 41688-0602 
283 99.9666 42.1221 $8.85641 -.24S0E-02 -. 3808E-02 -42778-02 
264 96.9666 60.2443 $8.46838 -.1708E-02 -. 3390E-02 - 3766E-02 
288 96.9994 $8 .3864 88.3790 -.1078E-02 -.1489£-03 . 1088E-02 
258 98.0000 ©.0000 68.9412 - 3479 E-02 -.7382E-08 . 34768-0602 
287 98.0000 8.8000 68.8264 -. 34328-0602 *. 7388 E-03 . 35810E-02 
. 288 96.0000 17.0000 58.8787 -.3362E-02 -.1462E-02 . 3986E-02 
289 66.0000 28.8000 86.8019 -.3114E-602 -.2182E-02 - 3788E-02 
260 66.0000 34.0000 86.6988 -.28118-02 -.2638EF-02 - 36968-02 
291 66.1661 36.0884 86.6336 - .26818-02 - .3087TE-02 - 4024E-02 
i 292 98.3122 42.1166 88.8665 -.2380E-02 - . 3347E-02 - 4090E-02 
233 98.4983 49.1783 88.4665 -.2182E-02 - . 33068-0602 - 3961E-02 
264 68.6244 §0.23358 88.4286 -. 21138-0602 - .32328-02 . 3891E-02 
266 68.7808 84.2622 66.3736 -.24328-02 -.18368-02 . SO4TE-02 
266 68.9365 88.3806 8S .3806 -.26918-02 *.11218-03 .2894E-02 
267 100.0000 ©.0000 $8 .9074 -.3280E-02 -.1218E-04 -32608-02 
268 100.0000 17.0000 85.8460 -.3188E-02 -.1428E-062 . 34918-0602 
298 100.0000 34.0000 88.8985 -.2686E-02 -.274&7E-02 . 38218-0602 
270 100.3117 42.1029 88.$438 - . 22708 -02 -.3281&-02 . 3988E-02 
271 100.6233 §0.2088 $8.4068 -.229868-02 -. 34868-0602 41328-0902 
272 100.9350 86.3088 68.3178 -.@1890E-02 -. 3I608E-03 -4204E-02 
273 102.0000 ©o0.90000 $s .8754 -.3126E-02 -.18 108-08 -3126E-02 
274 162.0000 &8.$000 88.8895 - . 3OSTE-02 -.71088-03 -3178E-02 
| 275 102.0000 17.0000 88.8150 - .3042E-02 -.13688-02 ~3347E-02 
276 102.0000 28.8000 538.7423 - .2808E-02 -.2027E-02 -3461E-02 
277 102.0000 34.0000 $6 .6427 *.2821E-02 -.2989E-02 - 3664E-02 
278 102.1866 38.0446 $6 .66456 -.2338E-02 -.2806E-02 - 3TIOE-O2 
279 102.3111 42.0890 $8.6216 - .218$3E-02 -.3189E&-02 - 36832E-02 
280 102.4967 488.1335 88.4840 -.2208E-02 *,3302E-02 - 3I9TOE-02 
281 102.6222 §$0.1780 $6.3846 -.2284E-02 - 37318 °02 -4389E-02 
282 102.7778 64.2226 686.3102 -. 4409E-02 -. 2736E-02 -5169E-02 
283 102.9333 88.2671 88.2670 ~.6914E-02 -.41368-03 -S$928E-02 
| 264 104.0000 ©.0000 88.8461 -.2931E-02 -.3140E-084@ -2631E-02 
c 268 104.9000 17.0000 88.7885 -.2880E-02 -,1362E-02 - 3186E-02 
288 104.0000 34.0000 88.6182 -.2378E-02 -,2S74B-02 - 3808E-02 
287 104.3222 42.0445 §8.5009 -.20818-02 - ,.3O073E-02 - 3TOOE-02 B-9 
104.9444 §0.0680 88.3808 -, Z24893E-02 - .3893E-02 438 7E-02 
104.9667 88.1335 $8.1814 - 13388-0601 -37TZIE-O2 -1389E-01 
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SEEPAGE FACE COMPUTATION RESULTS 


66 NE 8s 113 


TOTAL DISCHARGE 
113 NE « 196 


TOTAL DISCHARGE 
198 NE =: 215 


TDTAL OISCHARGE 
2158 NE =s 300 


TOTAL DIGCHARGE 
300 NE 317 
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TOTAL DISCHARGE 
604 NE =« 621 
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